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Abstract 

We consider polynomials orthogonal on [0, oo) with respect to Laguerre-type weights 
w(x) = x a e~Q( x \ where a > —1 and where Q denotes a polynomial with positive leading 
coefficient. The main purpose of this paper is to determine Planchercl-Rotach type asymp- 
totics in the entire complex plane for the orthonormal polynomials with respect to w, as well 
as asymptotics of the corresponding recurrence coefficients and of the leading coefficients 
of the orthonormal polynomials. As an application we will use these asymptotics to prove 
universality results in random matrix theory. 

We will prove our results by using the characterization of orthogonal polynomials via a 
2x2 matrix valued Riemann-Hilbert problem, due to Fokas, Its and Kitaev, together with 
an application of the Deift-Zhou steepest descent method to analyze the Riemann-Hilbert 
problem asymptotically. 

1 Introduction 

We consider Laguerre-type weights, 

w(x) = x a e~ Q{x \ forxe[0,oo), (1.1) 
where a > — 1, and where 

m 

Q(x) = Y,QkX k , q m >0, (1.2) 

denotes a polynomial of degree m with positive leading coefficient. In case Q(x) = x, the weight 
w is the classical Laguerre weight. Since all the moments of w exist we have a sequence of 
orthogonal polynomials. We use p n (x) = p n (x; w) = r y n x n + • • • with leading coefficient j n > 0, 
to denote the n-th degree orthonormal polynomial with respect to w, i.e. 

poo 

/ Pk{x)Pn(x)w(x)dx = 5 kjn , for k,n e N. 

J 

The orthonormal polynomials p n satisfy a three-term recurrence relation of the form, 

xp n (x) = b n p n+1 {x) + a n p n (x) + 6 n _ip n _i(x), n G N, (1.3) 
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and the coefficients a n and b n -\ are called the recurrence coefficients. 

The goal of this paper is to determine the asymptotic behavior (as n — * oo) of the recurrence 
coefficients a n ,b n -\ and of the leading coefficient 7„, and to determine Plancherel-Rotach type 
asymptotics for the orthonormal polynomials p n in the entire complex plane. In addition, we 
use these asymptotics to prove universality results for Laguerre-type unitary ensembles. 

The weights considered in this paper fall into the class of Freud weights on [0, oo). These 
weights are also closely related to symmetric Freud weights on R by a simple transformation so 
that our results can also be applied to Freud weights of the form |a:|Te - ^ a ) on R with 7 > —1. 
There is a vast literature on asymptotic questions with respect to Freud weights and we will 
now briefly summarize the results which are already known with respect to weights 

The leading order asymptotics of the recurrence coefficients associated with symmetric 
weights |x| 7 e - ^^ > on R, where Q is a polynomial with non-negative coefficients, have been 
derived by Damelin [I]. This result can then be applied to get the leading order asymptotics of 
the recurrence coefficients associated to the corresponding weights 

Full asymptotic description of the orthogonal polynomials and related quantities are known 
in the following special cases. First, the classical Laguerre case, i.e. Q(x) = x and a > —1, is 
described in |29| . Next, for a = ±1/2 and Q a polynomial, the asymptotics are known from 
results of Deift, Kriecherbauer, McLaughlin, Venakides and Zhou in 10 j. Indeed, the weights 
| x |±i/2 e -Q(a') Qn [0,00) are related to the symmetric weights e~® <yX ' on R, which belong to the 
class of weights considered in |1U| . See also the references therein for older results. Finally, in 
the special case that Q(x) = q monomial, one can verify that the weights x a e qmX 

on [0, 00) relate to the varying weights x a e~ nqrnX m on R, which fall into the class of weights 
considered by the author together with Arno Kuijlaars in 23 j. In that paper, the main goal was 
to prove universality of spectral correlations of the associated random matrix ensembles at the 
hard edge A = 0. Explicit formulae for orthogonal polynomials were not derived in [22] but can 
be deduced from their proof. 

One strong motivation for deriving the asymptotics of the orthogonal polynomials comes from 
random matrix theory. Indeed, the Laguerre-type weights (jl.lj) are related to random matrix 
ensembles used in the statistical description of mesoscopic quantum systems with disorder. 
Depending on the symmetries of the physical system one distinguishes ten different classes of 
ensembles, see for example |2"1 1181 l3*U] and references therein. Besides the three standard classes 
(orthogonal, unitary and symplectic) there are seven novel classes, which have been introduced 
by physicists during the past ten years. In the Gaussian case the joint probability distribution 
of the positive eigenvalues (Ei, . . . ,E n ) £ R™ can be written in the following form for all the 
seven novel symmetry classes, 

j- II \E]-Elf(f[Eje- E ?)dE 1 ...dE n , (1.4) 

11 l<j<k<n \£=1 J 

where (3 £ {1, 2, 4}, 7 £ N = {0, 1,2,...} and Z n is a normalization constant. We note that for 
these ensembles the eigenvalues come in pairs ■iiEj, i.e. for every positive eigenvalue Ej there 
corresponds a negavitive eigenvalue —Ej. In the spirit of the universality conjecture in random 
matrix theory we replace the term e~ E e in (|1.4|) by e~Q{ E D where Q denotes a polynomial of 
the form (|1.2j) . A simple transformation Xj = E? then leads to the following probability density 
function, 

1 n 1 

F n (X 1 ,...,X n ) = j- \ X j- X kfHxr 1] e- Q ^\ for (Ai,...,A n ) £ R+. 

^ n l<j<k<n 1=1 

It is a beautiful and important observation (see for example (SJI21]) in random matrix theory 
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that for f) = 2, 



F n (Xi, . . . , X n ) = — det(K n (Xj, \k))i<j,k<m 

where K n is constructed out of orthogonal polynomials with respect to the weight function 
w(x) = i^ 7 " 1 ^ - ^' 1 '. Moreover, all related statistical quantities of interest can be expressed in 
terms of K n . Using our results on the asymptotics of the orthogonal polynomials we will prove 
in this paper that for (3 = 2 the local eigenvalue statistics have universal behavior (as n — > oo) in 
different regimes of the spectrum. By universal we mean that the behavior is independent of Q 
as long as Q is chosen to be a polynomial. In work in progress, together with Deift, Gioev and 
Kriecherbauer fH], we will use the results of this paper, together with the approach developed in 
00 EI]) to prove universality for (3 = 1,4. 

In order to get our results, we are inspired by the papers of Deift et al. They considered 

orthogonal polynomials with respect to varying weights e~ nV ^ on R where V is real analytic 
and has enough increase at infinity [0], as well as with respect to non- varying weights e~®^ 
where Q is an even polynomial with positive leading coefficient jlUj . See also |Jj for an excellent 
exposition. As in these papers, we will characterize the orthogonal polynomials via the well- 
known 2x2 matrix valued Fokas-Its-Kitaev Riemann-Hilbert (RH) problem [I2|> and apply the 
Deift-Zhou steepest descent method, introduced in ^3] and further developed in ^HEUE]) ^° 
analyze this RH problem (and thus also the orthogonal polynomials) asymptotically. 

The main difference between the orthogonal polynomials with respect to the Laguerre-type 
weights 1)1- 1|) and the weights e~ nV ^ and e~^^ on R, lies in the behavior near the fixed 
endpoint 0. Near this endpoint we have to do a local analysis using Bessel functions. This 
is analogous to the local analysis near the fixed endpoints ±1 of the modified Jacobi weight 
considered in [2^, see also [T5| . 

The present paper is organized as follows. In the next section we state our results. In Section 
3 we state the RH problem for orthogonal polynomials and apply the Deift-Zhou steepest descent 
method to analyze the RH problem asymptotically. Afterwards, we use the asymptotics of the 
solution of this RH problem to determine the asymptotics of a n , 6 ra _i and 7„ (in Section 4), to 
determine the Plancherel-Rotach type asymptotics for p n (in Section 5), and to prove universality 
results in random matrix theory (in Section 6). 



2 Statement of results 

It is well-known in the theory of orthogonal polynomials that there are two ingredients needed 
for the asymptotic (as n — > oo) description. First, we need the sequence of Mashkar-Rahmanov- 
Saff (MRS) numbers, which we will denote by (3 n . For the weigths we consider in this paper, (3 n 
is uniquely determined for n sufficiently large by the equation 
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fin 



2W„ Q ' ix) \l—z dX = n - <"> 

see Proposition 13.41 From now on, we assume n to be sufficiently large such that the MRS 
numbers exist. The other ingredient we need, is the equilibrium measure [27] of [0, oo) in the 
presence of the rescaled field V n (x) = ^Q((3 n x). This measure is defined as the unique minimizer 
of the following minimization problem, 

inf log j-L-dn(s)dfi(t) + J V n (t)dfx(t)j , (2.2) 
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where the infimum is taken over all probability measures on [0, oo). We will show in Section \'A. 21 
below, that the minimizer is of the form, 



1 1 — x 

dfjL n (x) = ip n (x)dx = —d—^—h n (x)x(o,i]dx, (2.3) 

where X(o,i] denotes the indicator function of the set (0,1], and where h n is a polynomial of 
degree m — 1, 

m— 1 

h n (x) = h n , k x k , (2.4) 

k=0 

with real coefficients h n ^ : given by (|3.17|) . ()3.13|) and (|1.2I) . and which have an explicitly com- 
putable power series in ra -1 /" 1 , see (|3.19j) . 

2.1 Asymptotics of the recurrence coefficients a n ,6 n _i and the leading coeffi- 
cient 7 n 

To state our first result we need the quantity £ n = 2 J °° tp n {y) log ydy — V n (0). 

Theorem 2.1 The recurrence coefficients a n and 6 n _i in the three-term recurrence relation 
(11.31) of orthogonal polynomials with respect to the Laguerre-type weight ()1.1|) . have the following 
asymptotic behavior, 

S^ + 5^ + ° (1/ " 2) ' <* <2 ' 6) 
The leading coefficients 7 n of the orthonormal polynomials with respect to the Laguerre-type 
weight (jl.lj) have the following asymptotic behavior, 



7n = /?n (n+f + ^e-!"S/^ 



7T 

1 _ ( Aa2 ~ l 4. 12« 2 + 24a + ll _ ^(1) \ 1 2 
8^(0) ^ 24^(1) ^(l) 2 / « 1 7 ' 



, as n — » oo. (2.7) 



E'ac/i of the number sequences (3 n , t n , h n (0), h n (l) and h' n (l) as well as each of the error terms 
have an asymptotic expansion in powers of n~ l / m which can be calculated explicitly. 

Remark 2.2 For the convenience of the reader we summarize here the leading order behavior 
of the quantities appearing in (|2.5j) - (|2.7|) : 

oo m 

n = nV-^OV^ /3<°) = (\m qm A m )- l/m , A m = H 

3=0 j=l 3 

oo 

i n = Y e (j) n- j / m , = 4 log 2, 



j=0 



2m - 1 



/,„i(.) y,' 1 "" ' '"■ " 

i=o 

oo 

fr n (l) = Y h^n-M™, U 0) = 4m, 



3=0 



h' n (l) = V h^n~ j / m , A (0) = -m(m - 1), 

3=0 3 
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As 




-5 <5 1-6 1 1 + 6 



Figure 1: The different asymptotic regions for p n (/3 n z) in the upper half-plane. 



see Proposition 13.41 Remark 13.141 and equations ()3.18j) - (|3.22|) . 

Remark 2.3 The special case Q(x) = q m x m : In this case, it will be clear from Remarks 13.51 
13.111 and 13.151 that the values of n~ l / m (3 n , £ n , h n (0), h n (l) and h' n (l) do not depend on n and 
we have f3 n = n 1 /™^, l n = h n (0) = hP\ K(l) = and h' n {l) = hP\ Furthermore, 
the asymptotic expansions of the error terms in the theorem are given in powers of 1/n rather 
than in powers of n~ l l m . 



2.2 Plancherel-Rotach type asymptotics for the orthonormal polynomials p n 

In order to state the asymptotic behavior (as n — * cxd) of the rescaled orthonormal polynomials 
Pn(Pnz), with (3 n the MRS number, we need to introduce some more notation. Let 

1 (z-lW 2 

Mz) = 7T-- Tji—h n (z), for zGC\[0,l]. (2.8) 

Throughout this paper we always take the principal branch of the power. It is clear that 
ip n}+ (x) = —ip n -(x) = tp n (x) for x £ (0, 1), where ip n is the density of \i n . 

Further, we will show in Sections 13.61 and 13.71 below, that there exist biholomorphic maps f n 
and f n in a neighborhood of 1 and 0, respectively, satisfying 

2 A 1 

-fn(z) 3 / 2 = —win / ip n (s)ds, for \z — 1| small and z ^ (—cxd, 1], (2-9) 

3 Jz 

2f n (z) 1 ^ 2 = —mn \ ip n {s)ds, for \z\ small and z ^ [0, cxd). (2-10) 
Jo 

See QTGRt . (l3~?^ and (CTBHIl for an explicit expression of f n , and see (f3~TJTjt . (j3~T?rl and (l3~^4l 

for the explicit expression of f n . 

A last function which we will need is the conformal map ip from C \ [0, 1] onto the exterior 
of the unit circle, 

<p(z) = 2(z- 1/2) +2z 1 / 2 (z - l) 1 / 2 , for z G C\ [0,1]. (2.11) 



Because of the symmetry p n (z) = p n (z) we only need to present the asymptotics of p n {f3 n z) 
in the closed upper half-plane C+. We state the asymptotic formulae in the four closed regions 
As, Bs, Cs and Ds, depicted in Figure ^ 



Theorem 2.4 Let p n be the n-th degree orthonormal polynomial with respect to the Laguerre- 
type weight (|1.1|) . There exists 5q > such that for all 5 G (0, do) the p n have the following 
asymptotic behavior (as n — > oo). 
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(a) For z in the outside region Ag, 

p n (f3 n z) = e^*> 



Tlfin 2z 1 / i (z - 

(b) For z in the bulk region Bg, 



exp -mn / 4> n (s)ds (1 + C(l/n)). (2.12) 



z) 1 / 4 



cos f —(a + 1) arccos(2z — 1) — irn J ip n (s)ds — —j (1 + 0(l/n)) 
+ cos {^( a ~ ^) arccos (2z — 1) — 7rn ^ ip n (s)ds — — ^ 0(1/ 
(c) For z m i/ie Airy region Cg, 

Pn(PnZ) = GM"t e |Q(/3nZ) 



A» ^ 1/4 (^ - 1) 1/4 

cos ( i(a + 1) arccos(2z - 1)) /n(z) 1 / 4 Ai(/ n (z))(l + 0(l/n)) 



i sm 



J(a + 1) arccos(2z - 1)) f n {,z)~ x ^ Ai'(/ n (*))(l + 0(l/n)) 



with Ai the Airy function. 
(d) For z in the Bessel region Dg, 



fn{ 



with J a the J -Bessel function of order a, and with 

1 TTCt 

Ci(z) = -(a + 1) arccos(2z — 1) — . 



(2.13) 



(2.14) 



Pn zV 4 (l - 2)1/ 4 

(1 + 0(1/™)) 

+ cos &(*).£ (2(-/„(*)) 1/2 ) (l + 0(l/n))] , (2.15) 



(2.16) 



All the error terms are uniform for 5 in compact subsets of (Q,5q) and for z £ Xg with X £ 
{A,B,C,D}. The error terms have an explicit asymptotic expansion in n^ 1 /' rn . 

Remark 2.5 

1. The function arccosz which appears in (|2.13j) . 1)2. 14|) and (|2.16j) is defined as the inverse 
function of cosz : {0 < Re z < ir} — > C \ ((— oo, —1] U [1, oo)), and is given by, 



arccos z 



ds 



z (l-s)i/ 2 (l + s)i/ 2 ' 



for z G C \ ((— oo, — 1] U [l,oo)). 
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2. Some of the expressions in Theorem 12.41 are not well defined for all z£l, For example: 
z~2^ arccos(2z — 1), ip(z)^ a+1 \ (z — l) 1 / 4 , f* ip n (s)ds, etc. For these expressions we then 
take the limiting values as z is approached from the upper half-plane. 

3. The function ijj n (s)ds which appears in 1)2. 12j) and 1)2.13)) is explicitly computable. The 
result is, see Remark 13. 141 below. 

I ip n (s)ds = —H n (z)z 1/2 (1 - z) 1/2 - - arccos z 1/2 , for z e C+, (2.17) 

Jl 27T 7T 

where H n is a polynomial with real coefficients of degree m — 1, given by 1)3.48 1) . 

Remark 2.6 The above theorem is in agreement with |29i Theorem 8.22.8] where the Plancherel- 
Rotach type asymptotics for the classical Laguerre polynomials were stated. 

Remark 2.7 The special case Q(x) = q m x m : In this case /3 n = n 1 l m {^mq m A m )~ 1 l rr \ cf. 
R,emark )3.5) and since h n is n-independent it follows from ()2.8)) that ip n , and thus also L t/j n (s)ds, 
is n-independent. The latter integral is explicitly given by 1)2.17)1 where H n (z) = ^h(z) with 
h defined in Remark 13.151 Furthermore, the asymptotic expansions of the error terms in the 
theorem are given in powers of 1/n rather than in powers of n~ i ' m . 

2.3 Applications in random matrix theory 

Here, we consider matrices taken from the random matrix ensembles which induce the following 
probability distribution on the positive eigenvalues (Ex, ■ ■ ■ E n ) £ R", 

P^(E l ,...,E n )d n E = ±- J] \E 2 -E 2 \ 2 (flEje-^) ) d n E, (2.18) 



l<7<fc<™ V=l 



where we use d n E to denote dE\ . . . dE n . As noted in the introduction, after a transforma- 
tion Xj = E 2 the eigenvalue statistics can be expressed in terms of a scalar 2-point ker- 
nel K n constructed out of orthogonal polynomials with respect to the weight function 

W(X) = 12" 1 ' c , 

n-1 



K n (x,y) = V w (. x )V w (y)^2pk{x)pk(y) 



k=0 

= V w ( x )V w iy) • ( 2 - 19 ) 

7n x-y 

For example, consider the probability P n (LiE < x) that the smallest eigenvalue Le = mm j(Ej) 
is less than or equal to x G E+. Obviously we have P u (Le < x) = P n (L\ < x 2 ), with 
L\ = mhij(Aj) and Xj = E 2 . Since it is well known in random matrix theory, see e.g. 0, that 
Pn{L\ < x 2 ) = 1— det(I— fC n ), where K n is the trace class operator with integral kernel K n (x, y) 
acting on L 2 (0,x 2 ) and where det(J — JC n ) is a Fredholm determinant, we then obtain 

P n {L E <x) = l-det(I-IC n ). (2.20) 

For the case Q(x) = x it is known that the behavior of K n is described by the sine kernel in 
the bulk of the spectrum jWH^Hl, by the Airy kernel at the soft edge of the spectrum [Hfl l3Tj . 
and by the Bessel kernel at the hard edge of the spectrum )16) 132) . In the next theorem, we 
state that this behavior persists for any polynomial Q of the form 1)1. 2)1 . 

Theorem 2.8 Let w be the Laguerre-type weight l)l.l|) and let K n be the kernel ()2.19)) associated 
to w. Then the following holds. 
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(a) The bulk of the spectrum: With ip n {x) = ^ — ^rh n (x) the density of the equilibrium 
measure fi n , 



A, 



nip n (x) 



Pn [ X + 



ntp n (x) 



n^ n {x) 



sin(7r(ti — v)) 



+ 0[-\, (2.21) 



ir(u — v) 

as n —* oo, uniformly for u,v in compact subsets ofM. and x in compact subsets of (0, 1 
(b) The soft edge of the spectrum: With c n = (^/i n (l)) 2 ^ 3 , 



On.Tl 



2/3 



II 



c n n 



2/3 



Ai(u) Ai'(t;)-Ai(u)Ai'(«) 



+ o 



u — v 



nV3 



(2.22) 



as n — > oo, uniformly for u,v in compact subsets ofR 
(c) T/ie /tard edge of the spectrum: With c n = (|/i n (0)) 2 , 



4c n n 2 



4c n n z 4c n n z 

J a {y/u)^/vJ' a {^/v) - J a {^/v)^uJ' a {y/u) 



2{u-v) 

as n — > oo, uniformly for u, v in bounded subsets of (0, oo) 



n 



(2.23) 



Remark 2.9 As in |22| Section 3.4] one can show, using 1)2. 20|) and 1)2.23)) . that the smallest 
eigenvalue Le = min(£i, . . . , E n ), where the Ex,...,E n have probability distribution 1)2.18)1 . 
satisfies 



lim P n 



2cl/ 2 n 



ft 



1/2 



L E < s 



1-det [I- Ji 



;(7-l),^ 



2(«— y) 



ac- 



Here, we use J QjS to denote the integral operator with kernel 
ting on -L 2 (0, s). Tracy and Widom [32] have determined an explicit expression for the Fredholm 
determinant in the latter equation in terms of a Painleve V transcendent. Let q a (s) be the 
unique solution of the differential equation, 

~ l ){sq' a )' = q a (sq' a ) 2 + ~(s - a 2 )q a + ^sql{q 2 a - 2), ('= -^), 



with the boundary condition, 
1 



2 a r(l + a) 



a Q/2 



as s — > 0. 



Note that this equation is reducible to a special case of the Painleve V differential equation, see 
for example (22- Further, let 



F a (s) = exp f~ J log ^q a (x)dx 
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Then, it has been shown in [H2] that det(7 — J W)S ) = F a (s), so that 
2cl/ 2 n 



lim P n -^Le < s = 1 - Fi (7 _ 1} (^). (2.24) 



n— >oo 



Another application which we will present here concerns averages of products and ratios of 
characteristic polynomials det(xl — M) of random matrices M taken from the random matrix 
ensemble which induces the probability distribution (|2.18|) on the positive eigenvalues E\,...E n . 
Recall that the eigenvalues of M come in pairs ±Ej, so that det(xl - M) = ]lj=i0 2 ~ E j)- 
The averages are then intimately related with kernels Wi >n , Wii >n and Wm,n constructed out 
of orthogonal polynomials (with respect to the weight function w(x) = x^^^^e^^^') and 
their Cauchy transforms, see [3J I28j . Which kernel we have to use depends on whether the 
characteristic polynomials are only in the numerator (Wi jn ), only in the denominator (Wni,n) 
or in both the numerator and denominator (Wjx n ). We will not state for brevity reasons the 
explicit formulae of this connection here, but refer the reader to [HI I2H] (see also (S31 ) • The three 
kernels are given by, 

Wi n {u,v) = , (2.25) 

u — V 

WllnM = C7(7rnri>)(u)7r n . 1 (^-C(7r w - 1 u>)H7r B (i;) i ^ 

and 

C(7r n w)(n)C(7r n _iw)(-;;) - C(7r ra -m>)(M)C(7r n w)(u) 

Wni tn {u,v) = , (2.27) 

u — V 

with 7r n = ~fn l Pn = z n + ■ ■ ■ the ro-th degree monic orthogonal polynomial with respect to w, 
and C f the Cauchy transform of / on the positive real line, i.e. 

1 f°° f(r) 

Cf(z) = / ^-^dx, for z£C\ [0, oo). 

2-ki J x — z 

In the next theorem, we state that these kernels (and thus also the averages of characteristic 
polynomials) have universal behavior, as n — > oo, at the hard edge of the spectrum in terms of 
Bessel functions. 

Theorem 2.10 Let w be the Laguerre-type weight Ijl.lj) and let Wi^n,Wu n and Wm n be the 
kernels l|2.25jl - (|2.27|) associated to w. Further, let 3 a ,I> $alV ■ • ■ ^ e ^ e Vessel kernels given by 
Tabled and let c n = (^/in(0)) 2 . Then the following holds. 
(a) The kernel W\^ n satisfies, 

7n-l JrZ 9 W l,n [Pn ^r 2 '^"~ 



4c n n 2 ' \ 4c„n 2 Ac n n 2 

Pn \ ° O(0) ( tt t.. ,.\ , /n ( 1 



4c n n 2 

as n —* oo, uniformly for u,v in compact subsets of C 
(b) The kernel Wi\ n satisfies, 



e« >(j t ,,i(tM>)+o(i)), (2.28) 



as n — > oo, uniformly for u and v in compact subsets of C+ and C, respectively. The 
corresponding behavior of W\i^ n that holds uniformly for u and v in compact subsets of C_ 
and C, respectively, is stated by replacing with Jn a - 
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Limiting Bessel kernels 




„-«„-# Ja(uVa)«i/aj;(«Va)-j a ( w i/3)„i/aj;(„Va) 

2(u— d) 




Hi 1 ' ( U V2)„l/2j;(„l/2)_J Q („l/2 )li l/2 (if a) ) / (u l/2 ) 

a , 1 - 9 J^(t.V>i/V tt ( B VJ)-j a ( 1) i/>) tt i/>(4>))'( ! ,i/>) 
n2w 2 4(«-«) 


J a,IIl( M ' U ) 
J a,IIl( M ' U ) 


..*.j^ ) (« 1/ > 1/J (4 1, y(^ J )-fflf ) (» 1 / I )u l / a (fli 1) )'(« 1 / 8 ) 

U U 8(u-i;) 

u S u #i/«( U V2 )u l/2 {if W y(v l/ 2) _ ff P) (ul/2V/2{ff (l) ) , (ul/2) 

8(ii — vj 

fl a fl< 3 > («V> 1/3 (H (2) ) , ( „l/2 ( „l/2 )M l/2 (Jf (2) y^l/aj 

8(u— «) 



Table 1: Expressions for the limiting Bessel kernels. Here, J a is the usual J-Bessel function of 
order a, and Ha an d H^f 1 are the Hankel functions of order a of the first and second kind, 
respectively. 



(c) The kernel Wni.n satisfies, 

-(^&) % ^(^ , ' ) + G))■ (2 ' 30) 

as n — > oo, uniformly for u and v in compact subsets of C+ and C_, respectively. The 
corresponding behavior ofWm n that holds uniformly for u,v in compact subsets o/C+ is 
stated by replacing J^i a ^nia- Further, the behavior of Wni <n that holds uniformly 
for u, v in compact subsets of C_ is stated by replacing J in a with JTJJj Q . 

Remark 2.11 From our analysis we can also prove that the three kernels W\^ n , W\i n and Wm n 
(and thus also the averages of characteristic polynomials) have universal behavior at the soft 
edge 1 of the spectrum, and in the bulk of the spectrum. The proof of these facts is similar 
to the proof of Theorem 12.101 In the bulk of the spectrum the results will be analogous to the 
results of Strahov and Fyodorov in [2H] . 

2.4 The varying weights x a e~ nVt & 

Consider the varying weig hts x a e~ nV( - x) on [0, oo), where V is real analytic and has sufficient 
growth at infinity. We now briefly explain how one can obtain the analogues of the results of 
this paper for these weights. There are two ways to proceed. 

One approach is to link the orthogonal polynomials via a simple transformation to orthogonal 
polynomials with respect to symmetric weights of the form Ix^e - ™^^ ) on M. The latter weights 
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fall into the class of weights considered in |23| . Recall that in [2] the RH approach was used to 
prove universality for the associated random matrix ensembles, and not to determine asymptotics 
of the orthogonal polynomials. However, one can use the techniques of [221 to determine the 
asymptotics of the orthogonal polynomials. 

If one is interested in applications in random matrix theory, a second approach is more 
convenient since the kernel K n , see (|2.19[) . associated to the weig ht x a e~ nV ^ on [0, oo) doesn't 
relate directly to the kernel associated to the weight \x\^ e~ nVi - x > on R (they only relate via 
the orthogonal polynomials). One can apply the RH approach directly to the varying weights 
on the half real line. The analysis is then very similar as in the present paper, although the 
situation is more complicated. The complications come from the fact that the support Sy of 
the equilibrium measure \iy of [0, oo) in external field V can consist of more than one interval, 

L L 

Sy = [0, at] U (J [bj, a j+ i] or Sy = {J[bj,a j+1 }, a,,6j>0, 

3=1 3=1 

and that \xy can have degenerate behavior. 

One distinguishes four types of degenerate behavior depending at which point in R + the 
degenerate behavior occurs, cf. 0120]. Type I: points in [0, oo) \ Sy where equality holds in 
the variational condition (|3.1fi[) below. Type II: interior points of Sy where the density ijjy of 
fj,y vanishes. Type III: edge points ctj, bj where ij)y vanishes to higher order than a square root. 
And type IV: the edge point if E Sy and Vv(0) = 0. Near each of these singular points we 
have to construct special local parametrices. For type I-III singularities, the existence of these 
parametrices is established in U Section 5]. 

In the nondegenerate case, which is generic, the construction of the parametrices at the 
endpoints is as in the present paper. In particular, for regular soft edge points dj, bj this will be 
done using Airy functions, and for a regular edge point using Bessel functions. 

If Sy consists of more than one interval, the construction of the so-called parametrix for the 
outside region (see Section 1331 below for the one- interval case) uses ©-functions as in [SJ Section 
4.2]. 



3 Asymptotic analysis of the RH problem for orthogonal poly- 
nomials 

Here, we will perform the asymptotic analysis of the Fokas-Its-Kitaev RH problem for orthogonal 
polynomials ^3]. This RH problem is the following. Seek a 2 x 2 matrix valued function 
Y(z) = Y(z;n,w) which satisfies the following conditions. 

RH problem for Y: 

(a) Y :C\ [0, oo) -> C 2x2 is analytic. 

(b) Y possesses continuous boundary values for x £ (0, oo) denoted by Y + (x) and YL (x), where 
Y + {x) and Y-{x) denote the limiting values of Y(z) as z approaches x from above and 
below, respectively and 

Y + (x) = Y-(x)r X \ J, forx€(0,oo). (3.1) 

(c) Y has the following asymptotic behavior at infinity, 

Y ( z ) ( Z 7 °n) = 1 + 0(V*), as z ^ oo. (3.2) 
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(d) Y has the following behavior near z = 0, 



Y{z) 



O 



O 



o 



1 z a \ 

i \ 

1 log z 

1 log z ; 

: :)■ 



if a < 0, 



if a = 0, 



if a > 0, 



(3.3) 



as z -> 0, 2 G C \ [0, oo). 

Remark 3.1 The 0-terms in (j3,3jl are to be taken entrywise. So for example Y(z) = O (\ *Z) 
means that Y n = (9(1), Y 12 (z) = 0(z a ), etc. 

Remark 3.2 In contrast to the case of weights on the whole real line, considered by Deift et 
al. we now have an extra condition (d) near the origin, cf. ^HIE]. This condition is used 

to control the behavior near the origin. 



The unique solution of the RH problem for Y, see ^\ (for condition (d) see ^31^), is then 
given by, 



Y(z) 



TMz) 



hC(PnW)(z) 



In 



-27TZ7 n _ip n _i(z) -27rz7 n _iC(p n _iw)(2;) 



for zEC \ [0,oo), 



(3.4) 



where p n (x) = J n x n + • • • is the n-th degree orthonormal polynomial with respect to the weight 
w(x) = x a e~Q( x \ where 7„ > is the leading coefficient of p n , and where Cf(z) is the Cauchy 
transform of / on the positive real line, 



Am Jq s — z 



for z G C\ [0,oo). 



Remark 3.3 From (|3.4|1 . we immediately see that the kernels K n , Wi >n , Wu >n and Wni jn , given 
by ()2.19() and l|2.25|) - (|2.27|) . can be written in terms of Y. The kernels K n and Wi >n depend 
only on the first column of Y, the kernel Wn, n on both the first and the second column, and the 
kernel Wni,n only on the second column. Using the fact that detl" = 1, it is easy to check that 

K n (x,y) = x f e-^^yf e -5^)_J_-^ ( l) Y-\y)Y{x) Q , x,y € R, (3.5) 



and 



Wu tn (v,u) Wiii yn {u,v) 

-Wl, n (u,v) -Wll, n (u,v) 



-2irvy*_. 1 (u -v) 



Y- l (v)Y(u), u,v G C\ [0,oo). (3.6) 



Note in ()3.5() that (0 l) Y 1 and Y I J have an analytic continuation to 



Now, we will do the asymptotic analysis of the RH problem for Y. As in |3E3, see a ^ so !2]> 
we will use the Deift-Zhou steepest descent method |13| . and apply a series of transformations 
Y>-^Ut-^>Th^S^R to arrive at a RH problem for R with jump matrix uniformly close 
to the identity matrix. Then, one can show TO, that R is also uniformly close to the identity 



12 



matrix. By going back in the series of transformations we then have the asymptotic behavior of 
Y in all regions of the complex plane. 

The transformations are analogous to the corresponding ones in Yet, there are some 

technical differences which come from the factor x a in the weight. As noted in the introduction, 
the main difference lies in the fact that we have to construct a parametrix near the origin out 
of Bessel functions. The construction of this parametrix is analogous to the construction of the 
parametrix near the endpoints of the modified Jacobi weight, as done in |2j, see also |19| . 

3.1 MRS number (3 n and rescaling: Y \— > U 

The first step in the asymptotic analysis of the RH problem for Y will be a rescaling Y i— ► U. 
To do this rescaling, we will use the MRS number (3 n satisfying equation (|2.1|) above, and which 
will be constructed in the following proposition for sufficiently large n, cf. |1U1 Proposition 5.2]. 

Proposition 3.4 There is n\ G N such that for all n > n\ there exists a constant f3 n G K 
satisfying the condition, 



1 



2Wo Q ' {x) \j^—x dX = n - (3 ' 7) 
The number (3 n has a convergent power series of the form 

oo 

P n = n 1 / m Y,P [k) n~ k/m , (3.8) 
k=0 

with coefficients P^v that can be expressed explicitly in terms of the coefficients qo, . . . , q m of the 
polynomial Q. In particular, the first two coefficients [3^ and 0^' are, 

ft® = (lm qm A m y 1/m , /JW = - 2( 7" 1) ^- 1 , (3.9) 
2 ' m(2m — l)q m 

where 

A m = - f i m -\& = J] 2 ~^ 1 - (3-10) 
vr J \ l-x ajl 2j 

Proof. The construction of (3 n such that it satisfies ()3.7() is analogous to the construction of 
the MRS numbers a n ,(3 n in 10, Proposition 5.2]. Introduce the auxiliary function, 

1 f 1 m I t m 1 

G((3,e) = — ^^/^-V-M—dx = Y^^kq k A k f3 k e m - k . (3.11) 

n k=0 ' X k=0 

An easy calculation shows that, 

GiPn- 1 ^,^ 1 ^) = J-^Q'( x ) ] T^ldx. 

So, we need to construct (3 n such that G{[3 n n~ l / m , n _1//m ) = 1. 
Now, with /?(°) defined in (|33|) . we have by (|3~TT|) 



G(/3(°),0) = l, and ^(0,0) 



= m 

p=p<S>) 



^mq m A m ) 1/m / 0. 



Therefore, using the implicit function theorem, there exists £o > and a real analytic function 
p : (-e , e ) -» R, such that G(fi(e),e) = 1. We then define /? n = n 1 / m /?( n - 1 /™) for n sufficiently 
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large such that ra" 1 /™ < eq, and we obtain G(/3 n n _1//m , n -1//m ) = 1. We now have constructed 
Pn (for large enough n) to satisfy (|3.7jl . 

Since the function /3 is analytic near 0, it follows that n = n l / m [3( n - l / m ) has a convergent 
power series of the form (f3~%|) . The coefficients , /?( 2 ) , . . . can be expressed explicitly in terms 
of the coefficients qo, . . . , q m by an inductive argument using the fact that 



,, 2 

e-O fc=0 



for j = 1,2, 

e=0 



For example, for j = 1 this gives m 2 q m A m f3^ + (m — l)g m _i>l m _i = 0, so that [5^ is given by 

(EUD - □ 

Remark 3.5 The special case Q(x) = q m x m : In this case /3 n takes a simple form and exists 
for all n G N. By a straightforward calculation one can verify that f3 n = n 1 / m {^mq m A ra )^ 1 / m 
solves (|3.7j) in this special case. 

In order to rescale the RH problem for Y, we also have to introduce the rescaled field V n . 
Define, for all n > n\, 

1 m (\ \ 

V n (x) = -Q(PnX) = V -q k (3 k n X k . (3.12) 
n ^-^ \ n J 

k=o v 7 

So, ^ is again a polynomial of degree m with real coefficients. The coefficients of V n are 
n-dependent and have an explicitly computable power series in n -1 /" 1 . In particular one has 

m 1 oo 

V n (x) = Y^v n , k x\ v n , k = -q k fa = X: vfn { l m . (3.13) 

fc=0 l=m—k 

Remark 3.6 We immediately see, by using Proposition 13 A\ that the leading coefficient v nm 
has the asymptotic behavior v n ^ m = (^mA m ) _1 + 0(n _1//m ) as n — > oo, whereas the other 
coefficients of V n tend to zero as n — > oo. 

Remark 3.7 The special case Q(x) = q m x m : In this case V n (x) = ^q m {Pnx) m so that by 
Remark l3.5| y„(x) = (^m^m) -1 ^" 1 . Thus, w„ )m = (^mA m )~ l and u nj - = for j = 0, . . . , m— 1. 

Now, we are ready to rescale the RH problem for Y. Let (3 n be the MRS number constructed 
in Proposition 13.41 and define for all n > m, 

[/(*) = ^ (n+f),73 y(/3^)/3j a ' T3 , for z G C \ [0,oo), (3.14) 

with CT3 = ( o -i ) the third Pauli matrix. Then, it is straightforward to check, using (|3,12j) and 
the conditions of the RH problem for Y, that U is the unique solution of the following equivalent 
RH problem. 

RH problem for U: 

(a) U :C\ [0, oo) -> C 2x2 is analytic. 



(1 T a p -nV n (x)\ 
i )> forxG(0,oo). 

(c) £/(*) = (/ + 0(1/*)) f ^ , 
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(d) U satisfies the same behavior near z = as Y does, given by (j3,3j) . 



Remark 3.8 The RH problem for U is the RH problem for orthogonal polynomials correspon- 
ding to the rescaled weig ht x a e ~ nV ^) where V n is the rescaled field (l3~T2l) . 

The rescaling with (5 n is chosen to ensure that the equilibrium measure \x n of [0, oo) in the 
presence of the external field V n will be supported on the interval (0, 1). 

3.2 The equilibrium measure \i n of [0, oo) in the external field V n 

Here, we will determine the equilibrium measure fi n of [0, oo) in the presence of the external 
field V n . This measure will be used (via its log-transform) in the next subsection to normalize 
the RH problem for U at infinity. The equilibrium measure fj, n is the unique minimizer of (|2.2[) . 
and is characterized by the Euler-Lagrange variational conditions: there exists i G R such that 

2 J log \x — y\dfi(y) — V n {x) — 1 = 0, for x G supp( / u), (3.15) 

2 J log \x - y\dn(y) - V n (x) -£<0, for x G [0, oo). (3.16) 

The construction of \x n is as in jlUl Section 5.2]. It involves an analytic scalar function h n 
defined for all n> n\ as, 

K ™-hi.lj^M>Tr*' k!£C \ [al1 (3 ' 17) 

where V n is the rescaled field (|3.12|) . and T z is a positively oriented contour containing [0, 1] and 
z in its interior. This h n has the following properties, cf. .10, Proposition 5.3]. 

Proposition 3.9 The function h n is a polynomial of degree m — 1 with real coefficients that 
have an explicitly computable power series in n~ l l m . Furthermore, there exists n2 > n\ and a 
constant ho > such that h n (x) > ho for all n>n<i and x G [0, oo). 

Proof. From taking the residue at infinity in (|3.17|) we obtain after a straightforward calcula- 
tion that h n is a polynomial of degree m — 1 given by, 

m— 1 m 

h n(z) = K jk z k , h n , k = jvnjAj-k-l, (3-18) 

fc=0 j=k+l 

where the Aj are defined by (|3.1U|) and the v n j are the coefficients of the rescaled field V n . Recall 
that the v n j are real with an explicitly computable power series in n~ l l m . This yields that the 
h n ^ are also real with an explicitly computable power series in n -1 / m , 

j=U j=m—k 

Here, the leading order behavior of h n ^ has been determined by using ()3.18j) and Remark 
13.61 By (|3.19|) there exists n<i > n\ and a constant ho > such that for every n > i%2, all 
coefficients h n ^ are positive and h n (0) > Hq. Therefore, h n (x) > h n (0) > ho, for all x G [0, oo) 
and n > n<i- □ 
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Remark 3.10 Using (|3~T%j) and (|TT3|) we obtain that, 



m— 1 



A 



h n {z) = h(z) + 0{n- l ' m ), h(z) = 2 m ~ z k , (3.20) 



/, o Am 



as n — > oo, uniformly for z in compact subsets of C. Using 1, formula 15.4.1] the function h 
can be expressed in terms of a hypergeometric series, 

\vn 3 
h(z) = ^— [2 F 1 (l,-m + l;-m + -;z). (3.21) 

Furthermore, one can check (by induction on m) that 

4m 8 

h(0) = , h(l) = 4m, and h'(l) = -m(m - 1). (3.22) 

2m — 1 3 

Remark 3.11 The special case Q(x) = q m x m : In this case, it follows from ()3.18j) and Remark 
E2I that h n = h with h defined in ipOH]) . 

Now, we determine the equilibrium measure fj, n (for all n > 77-2) in terms of the polynomial 
h n , cf. ^3 Proposition 5.3]. 

Proposition 3.12 Define for all n > n 2 , 



1 1 — x 

*Pn(x) = —J —^—h n (x)x(o,i](x), (3.23) 

where h n is given by (|3,17|) . and where X(o,i] * s ^ e indicator function of the set (0,1]. By 
Proposition 13.91 it is clear that ip n is non-negative. Furthermore, 

1 

4> n (y)dy = 1, (3.24) 



and there exists a constant £ n S K smc/i f/ia£ 

2 y log |z - y\My)dy - V n (x) -£ n = 0, for x e [0, 1], (3.25) 
2 y log |z - y\My)dy - V n (x) -£ n <0, for x £ (1, 00). (3.26) 

So, dfi n (x) = tp n (x)dx is the equilibrium measure of [0,oo) in the external field V n . 

Proof. The proof is similar to [101 Proof of Proposition 5.3] and is based on the auxiliary 
scalar function 

Fn{z) = hi ^zS 2 hn{z) " 2^ y " (z) ' for z£C \ [0, 1]. (3.27) 

In order to prove that 4> n {y)dy = 1 we will determine two representations for the asymp- 
totic behavior of F n at infinity and compare them to each other. First, from the definition (|3.17l) 
of h n it follows that 

F n {z) = - iZ ~S /2 ±- f J^Kti)—, for . € C \ [0, 1], 
tti z v l l 2ir J \J I — y y — z 
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so that, since V^(y) = ^PnQ'iPnU) and by the condition IjM. 7j) on /3 n , 

1 1 /*• r~y~ r ,,, ,, , _ 2 , -1 



F -< 2 > - s; 2^ X V aT^' 3 '^ + 0(2 > = Si + 0(2 »• a8 2 - °°' (3 ' 28) 

Next, observe that F n<+ {y) — F n _{y) = 2tp n (y) for y G (0, 1). Since F n is analytic in C \ [0, 1] 
and since F n (z) = 0(z ) as z — » oo, see ()3.28|) . a standard complex analysis argument then 
shows that, 



A / X ^dy =±-.l l KAV) - FnAy) dy = F n ( Z ), forz€C\[0,l]. (3.29) 
m Jo y-z 2m J y-z 



Therefore, F n (z) = ^ J" Q tp n (y)dy + 0(z 2 ), as z ^ oo. Comparing this with (|3.28|) we obtain 
Jo 4>n{y)dy = 1. 

It now remains to prove the Euler-Lagrange variational conditions ()3.25[) and (|3.2(i|) . Since 
ip n G L 3 / 2 (M) it follows from and [H Theorems 5.31 and 5.32] that 

ImF n +(x) = — -f dy, a.e. for 

7T J x-y 

where the integral is a Cauchy principal value integral. Therefore, using also (|3.27|) . we obtain 
that 



d_ 

dx 



! J log \x - y\i> n {y)dy - V n {x)\ = 2-j^^^-dy - V^x) = 2vrImF n , + (x) - V^x) 

= - Re ({x - l)y 2 x; 1/2 /i n (x)) • (3-30) 

Since h n is positive on [0, oo), this yields that 2 J* log|x — y\ip n (y)dy — V n (x) is constant for 
x G [0, 1] and decreasing for x G (1, oo), so that ip n satisfies conditions (|3.25j) and (|3.26j) . □ 

3.3 Normalization of the RH problem at infinity: U i— > T 

In order to normalize the RH problem for U at infinity, we use the log-transform of the equi- 
librium measure d[x n (y) = if) n (y)dy. Define, for all n > n-i (with ri2 defined in Proposition 13.91 

above), 

9n(z) = / log (2 - y)i> n (y)dy, for z G C \ (-oo, 1], (3.31) 



o 



where we take the principal branch of the logarithm, so that g n is analytic in C \ (— oo, 1]. 

We now give properties of g n , cf. |1UI Proposition 5.4], which we will need in the following. 
From the definition of g n and from the Euler-Lagrange variational conditions (|3.25j) and ([3.26ft 
it follows that 

g n , + (x)+g n ,_{x)-V n {x)-£ n = 0, for xG [0,1], (3.32) 
2g n {x) - V n {x) - i n < 0, for x G (1, oo). (3.33) 

Furthermore, using ip n (y)dy = 1 we obtain 

9n,+ ( x ) ~ 9n-{x) = 2m, for x G (-oo,0), (3.34) 

9n,+ {x) - 571,-0*0 = 2vri / 4> n (y)dy, for x G [0, 1], (3.35) 

J X 
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and 



e n9 " iz) = z n + 0(z n - 1 ), asz^co. (3.36) 

Now, we are ready to perform the transformation U i— ► T. Define, for all n > 112, the matrix 
valued function T as, 

T{z) = e -h n ^U(z)e- n9 ^ z)e73 e^ n£na3 , for z G C \ R, (3.37) 

where £ n is the constant that appears in the Euler-Lagrange variational conditions (|3.25f) and 
(l3~26l) . Note that, by (l3~3^1 . the function e n9n has no jumps across (— oo,0), so that T has an 
analytic continuation to C \ [0, oo). It is then straightforward to check, using ()3.32j) . ()3.36|) and 
the conditions of the RH problem for U, that T is the unique solution of the following equivalent 
RH problem. 



RH problem for T: 

(a) T:C\ [0, oo) -» C 2x2 is analytic. 

(b) T + (x) = T^(x)vt{x) for x G (0, oo), with 
^ e -n(Pn,+(s)-ffn,-(a:)) 



vr(x) 







=«(ffn,+ (^)~3n,-(x)) 



I x a e n(2g n (x)-V„(x)-e n ) 

1 



for x G (0, 1), 
for x G [1, oo). 



(3.38) 



(c) T(z) = 1 + 0(l/z), as z -* oo. 

(d) T satisfies the same behavior near z = as F and J7 do, given by ()3.3j) . 

Remark 3.13 From (|3~35|) we see that the diagonal entries of vt on (0, 1) are rapidly oscillating 
for large n. From (|3.33f) the jump matrix vt on (1, oo) converges exponentially fast to the identity 
matrix as n — > oo. 



3.4 Opening of the lens: T ^ S 

Now, we will transform the oscillatory diagonal entries of the jump matrix vt on (0, 1) into 
exponentially decaying off-diagonal entries. This lies at the heart of the Deift-Zhou steepest 
descent method ^3], and this step is referred to as the opening of the lens. 

In order to perform the transformation T ^ S we will introduce scalar functions ip n and 
Define for all n > n-2, 

1 (~_ 1)1/2 

M?) = 77-^ T75— M*), for z G C \ [0, 1], (3.39) 

with principal branches of powers. So, the +boundary value of ip n on (0, 1) is precisely the 
density ip n of the equilibrium measure \i n . In particular we have, 

ip n ,+ (x) = -ipnA x ) = ^n{x), for X G (0, 1). (3.40) 
Now, define for all n > ri2, 

£ n (z) = -m J ip n (y)dy, for z G C \ (-oo, 1], (3.41) 
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where the path of integration does not cross the real axis. Note that in equation (5.34)] 
a function £ n is defined which is analytic through the support of the equilibrium measure. We 



find it more convenient to define ^ n with, ql branch cut along (0, 1). The definition of ^ n in (|3.41j) 
and ^3 equation (5.34)] also differs by a factor 2. 

The important feature of the function £ n is that, by (|3.40|) and (|3.35j) . £ ra)+ and £ n) _ are 
purely imaginary on (0, 1) and satisfy, 

2£„,+(x) = -2&, _(a;) = 2vri f My)dy = g n . + {x) - g n ,-{x), for x G (0, 1). (3.42) 

So, 2£ n and — 2£ n provide analytic extensions of g n ^ — g n - into the upper half-plane and lower 
half-plane, respectively. On R\ [0, 1), £ n satisfies 

in,+ {x) -£ n -(x) = 2iri, for x G (-oo,0), (3.43) 
2U(x) = 2g n {x) - V n (x) - £ n , for x G [1, oo). (3.44) 

Here, equation (|3,43|) follows from (|3.24|) . Equation (|3.44j) follows from the fact that the function 
2g n — Vn — 2£ n — £n is constant on [1, oo), by (|3.3U|) . and in x = 1, by (|3.25|) . Further, we can 
prove the existence of a 5± > such that for all n > ri2, cf. |10| Proposition 5.4] 

Re£ n {z) > 0, for < | Imz\ < Si and < Re z < 1. (3.45) 

We also need an estimate for £ n on (l,oo). From the fact that h n (x) > ho for all n > ri2 and 
x G [0, oo), see Proposition 13.91 we obtain 

1 ( x _ 1)3/2 

£ n (x) < --h -rjz — , for x G (l,oo) and n > n 2 . (3.46) 

o X ' 

Remark 3.14 As in jlUl Proof of Proposition 5.4 (v) and (vi)] we can determine £ n as well as 
the constant t n (which appears in the variational conditions) explicitly in terms of the coefficients 
v Uj k of V n . With arccos defined as an analytic function on C \ ((— oo, —1] U [1, oo)), as described 
in Remark 12 .5| one has 

£ n (z) =Ti(^H n (z)z 1/2 (l- z) 1/2 -2 arccos z 1/2 ^j , for z G C±. (3.47) 

where 

m— 1 / m \ 1 

R n(z) = E E VnjAj-k-i z k = -h(z) + 0{n~ l l m ), (3.48) 
k=0 \j=k+l J '" 

as n — > oo, uniformly for z in compact subsets of C. Here, the leading order behavior of H n has 
been determined by using Remark 13.61 and (|3.2()j) . 

Further, one can also verify that t n = — J2k"=o v n,k^k — 4 log 2. Then, since the coefficients 
v ni k have an explicitly computable power series in n~ 1 / m , so have the constants £ n , 

oo 

£ n = Y £^n- j / m , = 4 log 2, (3.49) 

j=0 

where the leading order behavior £^ of £ n has been determined by using Remark 13.61 

Remark 3.15 The special case Q{x) = q m x m : In this case, it follows from the previous remark 
together with Remark 13.71 that H n (z) = —h(z) and that £ n = —— — 41og2. 
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Figure 2: The lens shaped contour Ss = (Ji=i oriented from the left to the right. 



Inserting and (l3~4Hl into (jHSHJ, the jump matrix vt for T can be written in terms of 

the scalar function £ n as, 



v T (x) 



e -2ng nj _(x) 

1 ' 



for x G (0,1), 
for x G [1, 00) 



(3.50) 



A simple calculation, using the fact that £ n>+ (:r) + £ n _(x) = for x G (0, 1), see (|3.42|) . then 
shows that vt has on the interval (0, 1) the following factorization, 



v T (x) 



1 

x -a e -2n^ n -{x) 2 





-x~ 



X 



1 







, for x G (0,1), (3.51) 



and the opening of the lens is based on this factorization. 

Now, we are ready to do the transformation T 1— > S. Let Ss = u| =1 Sj be the oriented lens 
shaped contour as shown in Figure |^1 The precise form of the lens (in fact of the lips Si and 
S3) is not yet defined but for now we assume that it will be contained in the region where (|3.45fl 
holds. We will define the contour Sg explicitly in the beginning of Section f3.81 depending on n 
and on certain parameters 5 and v. Define, for all n > 712, an analytic matrix valued function S 
on C \ Sg as, 



S(z) 



(T(z) 



T(z) 



T{z) 



1 0" 



1 







for z outside the lens, 

for z in the upper part of the lens, 

for z in the lower part of the lens. 



(3.52) 



With the upper part of the lens we mean the region between Si and S2, and with the lower 
part of the lens the region between S2 and S3. 

One can easily check, using (|3.50|) . (J3.51|) and the conditions of the RH problem for T, that 
S satisfies the following RH problem. 

RH problem for S: 

(a) S : C \ S 5 -> C 2x2 is analytic. 
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(b) S+(z) = S-(z)vs(z) for z G Sg, with 



«s(z) 



1 1 

z -a g -2n^(2) j 

' z a \ 
-z- J' 

'l z a e 2n ^ n ^\ 
,0 1 ' 



(c) = / + 0(1/ z), asz^oo. 



for z G Si U S3, 
for z G S 2 = (0,1), 
for z G S4 = (1, 00). 



(3.53) 



Remark 3.16 Note that by ()3.45[) and H3.46j) the jump matrix vs on Si, S3 and S4 converges 
exponentially fast (as n — > 00) to the identity matrix. 

3.5 Parametrix p(°°) for the outside region 

From Remark 13 . 1 61 we expect that the leading order asymptotics of Y will be determined by a 
solution p(°°) , which will be refered to as the parametrix for the outside region, of the following 
RH problem. 



RH problem for p(°°): 

(a) P(°°) : C \ [0, 1] -> C 2x2 is analytic. 

(b) P^°\x) = p[°°\x) 



-x 



x° 

-a o 



for x G (0, 1). 



(c) P(°°)(z) =1 + 0(1/2 



as z — > cxd. 



As in |19| 1211 123j , we will construct a solution of this RH problem in terms of the Szego 
function D associated with x a on (0, 1). This is a scalar function which is analytic and non-zero 
in C\ [0, 1], that satisfies D + (x)D^(x) = x a for x G (0, 1), and which will not vanish at infinity. 
One can easily check that D is given by, 



D(z) 



y a/2 



if( Z ) a / 2 ' 



for z G C\ [0,1], 



(3.54) 



with principal branches of powers, and where <p is the conformal map from C \ [0, 1] onto the 
exterior of the unit circle, cf. (|2.11|) . 



ip(z) = 2(z - 1/2) + 2z 1 / 2 (z - l) 1 / 2 , for z G C \ [0, 1]. 
Since <p(z) = 4z + 0(1) as z — > 00, we have 
lim D(z) = 2~ a . 



(3.55) 



(3.56) 



The important feature of the Szego function is that the transformed matrix valued function 
2«o"3 _p(°o) d<t3 w qj satisfy conditions (a) and (c) of the RH problem, and that it will have the 
jump matrix ( _?i 0) on (0) !)• Then it is well known, see for example [H1E3]; that p(°°> is given 
by, 



p(oo)(. 



-00-3 



a(z)~ a(z) 1 a(z)—a(z) 1 

2 2~i 

a(z)—a(z)~ 1 a(z)+a(z)~ 1 

~i 2 



D(z)~ 



0-3 



for z G C\ [0,1], 



(3.57) 
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with 

(z - lW 4 

a(z) = ± — , for z G C \ [0, 1]. (3.58) 

Remark 3.17 Since p(°°) is a product of three matrices all with determinant one, we have that 
detp(°°) = 1. Further, for later reference, note that P ( -°°\z)z^ aa3 = 0(z -1 / 4 ) as z -> 0. 

Before we can do the final transformation S i— ► R we need to do a local analysis near 
and 1 since the jump matrices for S and p'°°) are not uniformly close to each other in the 
neighborhood of these points. 

3.6 Parametrix P n near the endpoint 1 

In this subsection, we will construct inside the disk Ug 2 = {zGC: \z — 1| < 62} with center 
1 and radius 82 > (sufficiently small and which will be determined as part of the problem in 
Prop osition 13 . 1 9l b elow) . a 2 x 2 matrix valued function P n that satisfies the following conditions. 

RH problem for P n : 

(a) P n : Ug 2 \ S 5 C 2x2 is analytic. 

(b) P„ j+ (z) = P n> _ (z)vs(z) for # G Sg n U§ 2 , with the jump matrix (|3.53|) for S. 

(c) P n {z)p(°°\z)- 1 =/+C(l/n) as n — » 00, uniformly for z on the boundary <9t/,5 of the disk 
Us and for 5 in compact subsets of (O,^)- 

The construction of P n has many similarities to the analogous construction carried out in 
[2lE]J) see also [Hj for an excellent exposition, and will be done using Airy functions. It typically 
involves three steps. First, we will construct a matrix valued function that satisfies conditions 
(a) and (b) of the RH problem for P n . In order to do this we will transform, in the first step, 

this RH problem into a RH problem for Pn^ with constant jump matrices and construct, in the 
second step, a solution of the latter RH problem. Afterwards, we will take in the third step also 
the matching condition (c) into account. 

Step 1: Transformation to constant jump matrices 

In order to transform to constant jump matrices, we seek the parametrix P n near 1 in the 
following form, 

P n {z) = E n (z)Pi 1 \z)e- n ^ z ^z- 1 2 a ^, for z G U &2 \ S s , (3.59) 

with E n an invertible analytic matrix valued function in U$ 2 , to be determined in the last step 
to ensure that the matching condition of the RH problem for P n is satisfied. 

The reader can easily verify, using ^ nj+ (x) +£„ ; _(x) = for x G (0, 1), that if Pn^ is analytic 
in Us 2 \ S5 with jump relations, 

for z G (Si U S 3 ) n Us 2 , 

y 1 1 1 

Pl%) = < 




Pi% )[ 2 1 ) • lor : € S 2 nr, 2 — ( \ — d : > . I ) . (:UiO) 

P nM z )y Q |V for zGS 4 nf/ 52 = (1,1 + 52), 

then P n defined by ()3.59j) satisfies conditions (a) and (b) of the RH problem for P n . In the next 
step we will determine P^ to satisfy these conditions. 
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Figure 3: The oriented contour 7 CT and the jump matrix v\ for $ on 7 ff . The four straight rays 
7(7,l> • • • ) 7cr,4 divide the complex plane into four regions I, II, III and IV. 



Step 2: Determine P„ explicitly 

The construction of Pn is based upon an auxiliary RH problem for \I> in the C-plane with jumps 
on the oriented contour j a , shown in Figure EJ consisting of four straight rays 

7<7,i : arg C = a, 7 CT ,2 : arg ( = vr, 7^3 : arg ( = -cr, 7^4 : arg ( = 0, 

with a E (? ,7r). These four rays divide the complex plane into four regions I, II, III and IV, also 
shown in Figure The RH problem for \& is the following, cf. HUj . 

RH problem for ^: 

(a) * : C \ 7 CT -> C 2x2 is analytic. 

(b) X f+(C) = (C) u i(C) f° r C £ 7(D where ui is the piecewise constant matrix valued function 
on 7 CT defined as shown in Figure 01 i.e. ui(C) = (J?) for £ G 7^1 and so on. This means 
that $ has the same jumps on 7 CT as pjp on Eg n Ug 2 , see (|3.60|) . 

(c) has the following asymptotic behavior at infinity, 



*(C) ~ C 



_Z2.J_ ( 1 1 



) k (sk-tk) Sk + tk 



_«0-o _2 C 3/2 



, (3.61) 



as C — > 00, uniformly for ( £ C \ 7 ff and a in compact subsets of (f ,7r). Here, 
r(3fc + l/2) . 6k + 1 



s k 



54 k k\T(k + 1/2) 



6fc - 1 



Sk, 



for fc > 1. 



(3.62) 
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It is well-known, see for example 0E3], that VP = VP CT (we suppress a in the notation for 
brevity) defined by, 

for C G I, 
for C G II, 

(3.63) 

for C G HI, 
for C G IV, 

with u> = e~ and Ai the Airy function, solves the RH problem for VP. See for example for 
definitions and properties of the Airy functions. 

Remark 3.18 Using 1, formulae 10.4.11 and 10.4.12] we have det VP = 1. This can also be seen 
from a different point of view as follows. From condition (b) of the RH problem for VP and from 
the fact that the Airy function remains bounded near it follows that det VP is entire. Further, 
from condition (c) we have det VP — > 1 as z — > oo. Using Liouville's theorem we then indeed 
obtain det VP = 1. 

The idea is now to construct pjp out of VP as Pn\z) = \P(/ n (z)) for appropriate biholo- 
morphic maps f n : Us 2 — > f n {U$ 2 ) with / n (l) = 0. We will choose these biholomorphic maps to 
compensate for the factor e~ ri £ n W <73 in Q3.59JI . So, by the asymptotic behavior Q3.61JI of VP at 
infinity, we need, 

-\fn(zf' 2 = nUz), for zeU d2 \ (-oo, 1]. (3.64) 

Note that, by (|3.41|) . this is precisely equation (|2,9|) . The construction of these biholomorphic 
maps is analogous as in |101 Section 7.1] and we define, for all n > ri2, 

f n {z) = n 2 / 3 n (z), for z G Ug 2 , (3.65) 

with (f) n defined in the following proposition, cf. 10, Proposition 7.3]. 

Proposition 3.19 There exists 82 > such that for all n > n,2 there are biholomorphic maps 
4>n ■ Us 2 -> 4> n {Us 2 ) satisfying: 

1. There exists a constant cq > such that for all z E U$ 2 and all n > ri2 the derivative of 
(j) n can be estimated by: Co < |0^(z)| < 1/cq and | arg0' n (z)| < 

2. 4> n {Us 2 n M) = MUs 2 ) n R, and 4> n (u S2 n C±) = MUs 2 ) n C±. 

5. -lMz) 3/2 = Cn(z) for zGU S2 \ (-00, 1]. 
Proof. Define, for all n > ri2, the auxiliary function, 

Mz)= 2 7^3/2 , for zGC\ (-oo,l]. (3.66) 



VP(C) = v^e § x < 



/Ai(C) Ai(u; 2 C) 

^Ai'(C) u; 2 AiVC), 

/Ai(C) Ai(u; 2 C) N 

{m>(() w 2 AiV0, 

/Ai(C) -^ 2 Ai(u;C) N 

\m'(0 -Ai'K), 

/Ai(C) -^ 2 AiK) N 

Ui'(C) -Ai'(wC) 



e e 



e e 



<T3 



cr.-j 



e e 



1 

; ::) 
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Note that by l|3.42j) the function cf) n has no jumps across (0,1), so that 4> n has an analytic 
continuation toC\((-oo,0]U{l}). From (J3~4T|) and it follows that, 



(3.67) 



Using Cauchy's theorem and the fact that h n is uniformly bounded in compact subsets of C, 
which follows from equation (|3.20j) . there exists a constant c > such that for all n > n2 and 
\s- 1| < 1/4, 



h n (s) 



A/2 



h n {l) 



1) 



1 

'27ri 



h n (w)w 1 I 2 - h n (l) dw 



\w-l\=$ 



W 



< 4|s — 1| sup 



h n (w) 



w 



1/2 



hjr 



i 

< els 



w 



II. 



|f-iN 

Inserting this into (|3.67j) we obtain that there exists a constant Ci > such that 
\4> n {z) — 1| < C\\z — 1|, for all n > n-2 and |z — 1| < 1/4. 



(3.68) 



Therefore, the isolated singularity of 4> n at 1 is removable so that <fr n is analytic in C \ (— oo,0], 
and there exists 5 > such that ~Re<f) n (z) > 0, for all n > ri2 and \z — 1| < 5. This yields, 



4>n{z)^{\h n {l)) 2/ \z-l)Uz? IZ 



(3.69) 



is analytic for z £ Us- 

Observe that, by (|3.68j) and (|3.69|) . 4> n {z) is uniformly (in n and z) bounded in Us- This 
implies, by using Cauchy's theorem for derivatives, that 4>n( z ) is also uniformly (in n and z) 
bounded in Us for a smaller 5. Since (f> n (1) = 1, see fLEEl, we have (f>' n (l) = (i/i n (l)) 2 / 3 , so that 



" (o^n(l)) 



2/3 



< C 2 \z-1\, 



for all n > n2 and z E C/^, 



for some constant C2 > 0. Therefore, since h n (l) > ho > 0, see Proposition 13.91 there exists 
< 62 < 8 such that for all n > ri2 the <j) n are injective and hence biholomorphic in U$ 2 and 
such that they satisfy part 1 of the proposition. 

Part 2 follows from the first part (for a possible smaller 62)- The last part of the proposition 
follows from the second part and from equations (|3.69l) and (|3.66|) . □ 

Remark 3.20 For later reference, observe that by (|3.65fl and (|3.69|) the biholomorphic maps 
/„ are given by, 

fn(z) = c n n 2/3 (z - l)fn(z), for z G U$ 2 , 



(3.70) 



where f n = (f) n with cf) n given by ()3.66|) . and where c n = (^h n (l)) 2 ^ 3 . The constant c n has, by 
Remark 13,101 the following asymptotic behavior, 

c n = {2m) 2 /\l + 0(n- l l m )), 



as n 



00. 



(3.71) 



Furthermore, from the proof of the proposition it follows that f n is analytic and uniformly (in 
n and z) bounded in Us, for some 5 > 82, and that f n iX) = 1- Therefore, there exists a constant 
C > such that 



- 11 = 

for all n> 17,2 and z E c7,5 2 



f n (s) - I ds 



\s-l\ = ^ 



1 



1| < C|z- 1|, 



(3.72) 
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We now have introduced the necessary ingredients to define Pn . Let n > ri2 and a G (f , vr), 
and recall that the contour £5 is not yet defined. We suppose that S5 is defined in Us 2 as the 
inverse / n -image of 7 CT n f n (Us 2 ). Define, 

P^(z) = *(/„(*)), for z G Ug 2 \ f~\ 1(T ). (3.73) 

Then, we immediately see that P n is analytic in Us 2 \ with jump relations (|3.60|) . 

Step 3: Determine E n explicitly 

In this final step, we determine the invertible analytic matrix valued function E n in equation 
(|3.59j) such that the matching condition (c) of the RH problem for P n is satisfied. From (|3.59|) . 
(|3.73|) . ()3.61j) . ()3.64j) and (|3.65|) we see that (to ensure that the matching condition is satisfied) 
we have to define E n , for all n > ni as, 

E n {z)=P^\z)z^e^ -L (] ~f\Mzft, foizeU S2 . (3.74) 

This ends the contruction of the parametrix P n . 

Remark 3.21 Obviously, E n is analytic in U$ 2 \{— 00, 1]. Using condition (b) of the RH problem 
for P(°°) and using the fact that / nj+ (x) <T3 / 4 = f n _(x) <T3 / 4 e^ CT3 for x G (1 - 8 2 , 1), it is easy to 
check that E n has no jumps on (1 — 82, 1). So, what remains is a possible isolated singularity at 
1. However, E n has at most 1/2-root singularities at 1, which implies that the singularity at 1 
has to be removable. Therefore, E n is indeed analytic in U§ 2 . 

Further, from (|3.74|) and from the fact that det-P^ 00 -* = 1, see Remark 13.171 it follows that 
det E n = 1, so that E n is also invertible. 

Summary of the obtained result 

We will now briefly summarize the obtained result. Let n > ri2 and a G (f ,tt), and suppose 
that the contour E5 satisfies f n (J^s H Us 2 ) = 7 CT H f n (Us 2 ). Define, 

P n (z) = E n (z)^(f n ( Z ))e-<^ Z -^, for z G Ug 2 \ f' 1 ^), (3.75) 

where the matrix valued function E n is given by (|3.74j) . the matrix valued function by 1)3.63(1 . 
and the scalar function f n by ()3.65(l . Then, P n solves the RH problem for P n . Furthermore, 
using (jTTBl . (jTTH . (JMTJ), (J3HJ), ^U^i and part 3 of Proposition EH we have 

00 1 

Pn^P^iz)- 1 ~J + ^A fc (z)-^, asn^cx), (3.76) 

k=l 

uniformly for z in compact subsets of {0 < \z — 1| < 62} and a in compact subsets of (§,7r), 
where A/% is a meromorphic 2x2 matrix valued function given by, 

x( ^ l ti? + t \\ l{Sk -' k) )z-^P^\zr\ (3.77) 
\-i(-l) k (s k -t k ) s k + t k J 1 ; v ; 

for z £ {0 < \z — 1\ < 82} ■ Here, the coefficients s k and t k are defined by (|3.62j) . 

Remark 3.22 Obviously, A k is analytic in Us 2 \(— 00, 1]. Using the fact that £„ )+ (x) = — £ n) _(a;) 
for x G (1 — 82, 1) together with condition (b) of the RH problem for p(°°\ the reader can verify 
that Afc has no jumps on (1 — 82, 1), so that A^ is indeed meromorphic in U$ 2 . 

Furthermore, since £ n (z) = 0{(z - l) 3 / 2 ) and P^°°\z) = 0({z - l) 1 / 4 ) as z 1 it follows 
that Afc has a pole of order at most [ 3fc 2 1 ' 1 ] at 1. 
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3.7 Parametrix P n near the endpoint 

Here, we do the local analysis near 0. We will construct inside the disk Us 3 = {z £ C : \z\ < 5%} 
with center and radius £3 > (sufficiently small and which will be determined as part of the 
problem in Proposition 13. 25l below). a 2 x 2 matrix valued function P n that satisfies the following 
conditions. 

RH problem for P n : 

(a) P n : Us 3 \ S 5 C 2x2 is analytic. 

(b) P nj+ (z) = P n -(z)vs(z) for z £ n Ug 3 , with vs the jump matrix ()3.53|) for S. 

(c) P n {z)P^(z)- 1 =I + 0(l/n), as n — ► co, uniformly for z E dU$ and for 5 in compact 
subsets of (0, £3). 

The construction of P n is similar to the construction of the parametrix near the endpoints ±1 
of the modified Jacobi weight, see |191 121j . and will be done using Bessel function. It involves, 
like the construction of P n , three steps. In the first two steps we construct a matrix valued 
function that satisfies conditions (a) and (b) of the RH problem for P n . In the last step we take 
also the matching condition (c) into account. 

Step 1: Transformation to constant jump matrices 

Seek P n in the form, 



with E n an invertible analytic matrix valued function in U$ 3 , which will be determined in the 
third step. Note that by ()3.43|) the function e~ n ^ n has no jumps across (— 00, 0). Therefore, if 
PrP is analytic in Us 3 \ so is P n . 

It is straightforward to check, using £n,+ 0*0 + Cn-(x) = for x G (0,1), that if P{p is 
analytic in Ug 3 \ S5 with jump relations, 



then P n defined by (|3.78|) satisfies conditions (a) and (b) of the RH problem for P n . 
Step 2: Determine PrP explicitly 

The construction of P^ is based upon an auxiliary RH problem for ^ a in the £-plane with 
jumps on the oriented contour j a , shown in Figure |1J consisting of three straight rays 

% tl : arg ( = -a, % >2 : arg £ = 7r, % >3 : arg £ = a, 

with a £ (0, 7r). These three rays are oriented to infinity and divide the complex plane into three 
regions I', II' and III', also shown in Figure 2] This auxiliary RH problem has been used before 
|19( 121 j in the construction of the parametrix near the endpoints ±1 of the modified Jacobi 
weight, and is the following. 



P n (z) = E n {z)P£Xz)e- n ^ z ^{-z)-\ 



for z e U Ss \ £5, 



(3.78) 




(3.79) 
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Figure 4: The oriented contour 7^ and the jump matrix vq for ^ a on j a . The three straight 
rays 7 CT) 1,7^,2 and 7 CT3 divide the complex plane into three regions I', II' and III'. 



RH problem for 

(a) ^ a : C \ % -> C 2x2 is analytic. 

(b) ^^^(C) = \& a) _(£)wo(C) for C G 7(j, where vo is the piecewise constant matrix valued 
function on % defined as shown in Figured i.e. i>o(C) = ( e -^<* 1) f° r C G 7 CT ,i, and so on. 

(c) has the following asymptotic behavior at infinity, 



✓ 4fc,Afc/2 



2 2?1/2<T3 , (3.80) 



as C — > 00, uniformly for ( £ C \ > and cr in compact subsets of (0,7r). Here, (a, 0) = 1 
and 



(a, k) 



(4a 2 -l)(4a 2 



(4a 2 - (2k - If 



2 2k kl 



for k > 1. 



(3.81) 



One knows |19| I21j that the 2x2 matrix valued function ty a = vl/^ (we suppress a in the 
notation for brevity) defined by, 



*«(C) 



I a (2(s) -^a(2C 5 ) 
-2^4(2^) -2(^(2^ 

i^(2(-C)^) 



<5( J ff( 1 ))'( 2 (-C)5) vrC^^ , )'(2(-C)5) 



|^i 2) (2(-C) 

(2). 



Q7rZO"3 



i^ 2) (2(-C) 



|^ 1} (2(-C) 



k Uc^(^ 2) )'(2(-C)^) <5( jF /( 1 ))'(2(-C)^ 



for C G I', 



for C E II', (3.82) 



e-3 aBri0 *, for C GUI', 
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solves the RH problem for $> a . Here I a and K a are modified Bessel functions of order a, and 
H ( a ] and H {2) are Hankel functions of order a of the first and the second kind, respectively. See 
for example for definitions and properties of these functions. 

Remark 3.23 From 1, formulae 9.1.17 and 9.6.15] we have det^ = 1. 

As in the construction of the parametrix near 1 we construct Pn out of ty a as ^ a (f n (z)) 
using appropriate biholomorphic maps f n : Ug 3 — > f n (U$ 3 ) with / n (0) = 0. We will choose them 
to compensate for the factor e~ n ^ z ' <T3 in (|3.78j) . By (|3.80j) we see that a good choice would be 
to construct f n such that it satisfies, 

e vT(z) = (_ 1 )n e ni„(*) j for z G U$ 3 \ [0, oo). (3.83) 

In order to construct f n we prove the following lemma and proposition. 

Lemma 3.24 Define, for all n > ri2, the auxiliary scalar function, 



1/2 l-7^ I Ms)ds, (3.84) 



-z 



h n (0) Jo 



with ip n given by (|3.39|) , This function is, by (|3.4U|) . well-defined and analytic in C \ [l,oo). 
Furthermore, there exist constants C\ , 63 > such that 

4> n (z) — 1 < Ci|z|, for all n > n2 and \z\ < 1/4, (3.85) 

Re ((-z) 1/2 4> n (z^J > 0, for all n>n 2 and z G \ [0, 00). (3.86) 

Proof. Equation (l3~85ll can be proven analogously as equation (|3.68|) in the proof of Proposi- 
tion cnni 

By ([3.85ft there exists 5 > such that | arg(fi n (z)\ < vr/4, for all n > ri2 and \z\ < 5. Since 
I arg(— z) 1 / 2 ) < 7r/4 for | argz| > ir/2, we then have 



Re \ {-z) 1/2 ^ n {z)J > 0, for all n > n 2 and \z\ < S with | argz| > f . (3.87) 
Now, from (|3^U|) . (|3^T|) and the fact that / ifi n {s)ds = 1, we obtain that 

, - ~ri fz I u \n\ (£n(z) — TTl), iflmz>0, 

(-z)^k(z) = / M')ds = *f ' ] (3.88) 

hn ^ Jo [li^^n(z)+m), iflmz<0. 

So, by (|3.45j) and since h n (0) > 0, 

Re ((-z) 1/2 n (z)) > 0, for all n> n 2 and \z\ < 8 X with < |arg^| < §. (3.89) 

with 5\ defined in (|3.45j) . From (|3.87j) and (|3.89|) . equation ()3.86jl is then proven with £3 = 
min{5, 5\}. □ 

Proposition 3.25 There exists 63 > such that for every n > ri2 there are biholomorphic maps 
4>n ■ Us 3 -> 4> n {U& s ) satisfying, 

1. There exists a constant cq > such that for all z G U~s 3 and all n > 77-2 the derivative of 
(fi n can be estimated by: cq < 1^0*01 < V c o an d I arg^nl 2 ) — 7r| < 7r/15. 
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2. e 2n4>i /2 (z) _ (.Jjn^W ^ or ze fj 6a \ [0, oo). 

Proof. Let 

4> n (z) = ~h n (0) 2 z4>l(z) = Q/ ln (0)(-z) 1 / 2 ^ n (z)^ , for z G C \ [1, oo), (3.90) 

with cfi n defined by (|3.84|) . From (|3.85|) and from the fact that h n (0) is bounded it follows that 
(fin{z) is uniformly (in n and z) bounded in {|z| < 1/4}. So, by using Cauchy's theorem for 
derivatives, <fi n {z) is uniformly (in n and z) bounded in {\z\ < 5} for some < 5 < 1/4. Since 

4> n (0) = 1, see (UTEHJ), we have 0' n (O) = -\h n (ti) 2 , and thus 



+ ^ n (0) 2 



< C2I2I, for all n > ri2 and |z| < 5. 



Therefore, since h n (0) > ho > 0, there exists < 63 < 5% such that for all n > n2 the (/> n 
are injective and hence biholomorphic in U$ 3 and such that they satisfy the first part of the 
proposition. 

The second part of the proposition can be verified by using equations ()3.86|) . <|3-88|) and 
(l3~90l . □ 

We now define the biholomorphic maps f n for all n > ri2 as, 

f n (z) = n 2 4> n (z), iorzeUs,. (3.91) 

By the second part of the proposition, equation (|3.83j) is then satisfied. Note that, by (|3.84|) . 
(ET%51) and (jSSOl), equation (EHUj) is satisfied. 

Remark 3.26 For later reference, we state the analogue of Remark 13.201 From (|3.90p and 
(GSHJ), it follows that 

f n (z) = -c n n 2 zf n (z), for z G Ug 3 , (3.92) 

where f n = 4> n with cj) n given by ()3.84() . and where c n = (^h n (0)) 2 . The constant c n has, by 
Remark 13.101 the following asymptotic behavior, 



r, \ 2 

2m 



2m -1 



(1 + 0(n- 1 ' m )), asn^oo. (3.93) 



Furthermore, from Lemma 13.241 we have that f n is analytic in Ug 3 and that / n (0) = 1. As in 
Remark 13.201 there exists a constant C > such that 

fn(z) - 1 <C\z\, for all n > n 2 and z G Ug 3 . (3.94) 

Remark 3.27 Observe that, in contrast to the biholomorphic maps (fi n of Proposition 191 the 
function cfi n maps the upper (lower) part of the disk Ug 3 onto the lower (upper) part of (fi n (Us 3 ). 

We now have all the ingredients to define . Let n > n 2 and <r G (0,7r), and recall that 
the contour S5 is not yet defined. We suppose that E5 is defined in Us 3 as the inverse / n -image 
of jo- n f n (Us 3 )- Define, 

Pjp{z) = * a (j n {z)), for z G Us 3 \ fn\%)- (3-95) 

Then, is analytic in Ug 3 \ S5 with jump relations l|3.79|l . 
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Step 3: Determine E n explicitly 

In this final step, we determine the invertible analytic matrix valued function E n in equation 
l)3.78j) such that the matching condition of the RH problem for P n is satisfied. From ()3.78j) . 
(I3~95l) . iprSH]) . (l3~83l) and (l3~9T]l we have to define E n , for all n > n 2 as, 

E n (,) = (-irpH(,)(-,)>a^l j)/ B (*)^(27rr/ 2 , farzeET*. (3.96) 

This ends the construction of the parametrix near 0. 

Remark 3.28 Obviously, E n is analytic in Ug 3 \ [0, oo). From condition (b) of the RH problem 
for p(°°) and from the fact that f n ^(x) a3 ^ 4 = e~~ a3 f n -{x) a ' i ^ for x £ (O,^), one obtains 
that E n has no jumps on (0,53). So, what remains is a possible isolated singularity in 0. 
However, using Remark I3.17( one sees that we have at most 1/2-root singularities in 0, so that 
the singularity at the origin has to be removable. Therefore, E n is indeed analytic in Ufa- 

From (|3.96|) and from the fact that det p(°°) = 1, see Remark 13.171 it follows that det E n = 1, 
so that E n is also invertible. 



Summary of the obtained result 

We will now briefly summarize the obtained result. Let n > n 2 and a £ (0, tt), and suppose that 
the contour S5 satisfies / n (Ss fl Ug 3 ) fl fn(Ug 3 ). Define, 

P n (z) = E n (z)^ a {f n {z))e~<^(-z)-^\ for z£Ug 3 \ J-\%), (3.97) 

where the matrix valued function E n is given by ([3.96ft , the matrix valued function ^ a by (|3.82|) , 
and the scalar function f n by (|3.91() . Then, P n solves the RH problem for P n . Furthermore, 
using (JinZZJ), d3~96j) . (l3~80|) . (l3~83|) and (|3~9TTl . we have 

00 1 

P^P^z)- 1 ~/ + VA,.(zW, asn^cx), (3.98) 

^-^ n K 

k=l 

uniformly for z in compact subsets of {0 < \z\ < 63} and for a in compact subsets of (0, tt), 
where A/% is a meromorpic 2x2 matrix valued function given by, 

x ( ^(« 2 + (k - \)i \ i aCT3p(oo) x 

for z £ {0 < \z\ < 63}. Here, (a, 0) = 1 and (a,k) is defined by (|3.81j) for k > 1. The function 
4> n is defined by (j33cTl) and 

Remark 3.29 We have the analogue of Remark l3.221 Namely, one can check that A^ is indeed 
meromorphic in Ug 3 and has a pole of order at most [^jr 1 ] at 0. 

3.8 Final transformation: S 1— > R 

In this subsection, we will perform the final transformation of our RH problem. Recall that the 
contour Eg is still not yet explicitly defined. We will now define it in terms of the parameters 
n, 5 and v (a new parameter replacing a). Here, we follow Section 7.2]. 



31 




Figure 5: The contour £# = Eg U dlf$ U dU$ depending on the parameters n, 5 and v. Here 
£ 5 = Ui=i s i with s i = s 'i u s i u s i' and so on- 



Let 5 = min{(5i,52,53} (cf. (j23SJ), Proposition l3~TT)l and l3~25)) . Fix 5 £ (0,5 ), n > n 2 , 
(cf. Proposition 13.9(1 and v 6 From Proposition we know that there exists a 



cr(n, z/, 5) G ( 



2_77 

3 



15' 6 



+ such that /-Htt.i) n dU s 



{1 + 5e tu }. By the symmetry 
/nO 2 ) = /n(^) we then also have f^ila^) H 3C/<5 = {1 + 5e~ lu }. We then define Eg in C/5 as 
the inverse / n -image of ^/o-. We can do an analogous construction near 0, and define Eg in the 
disk Us near as the inverse / n -image of ^a(n,v,8) such that \la,l) D dU$ = {— 5e~ lu } and 

f- 1 (% 3 )ndu s = {-Se iv }. 

Further, define a contour £# in terms of the contour Eg. Let £# = Eg U df/^ U dU$. This 
leads to Figure |SJ Note that the contour £# depends on n (and also on 5 and z^). However, we 
immediately see that £'", . . . , £4' are independent of n. 

Now, we are ready to do the transformation S *—> R. Define a matrix valued function 
R : C \ Sij ; — » C 2x2 (depending on the parameters n, 5 and v) as, 



S^P^)- 1 , 
S(z)P n (z)-\ 
S(z)P(°°\z)-\ 



for z £ Us \ Eg, 
for z £ Ug \ Eg, 
for z elsewhere, 



(3.100) 



where P n is the parametrix near 1 given by (|3.75j) . P n is the parametrix near 0, see 

is the parametrix for the outside region given by (|3.57() . and S solves the RH problem for S. 

Remark 3.30 The inverses of the matrices P n , P n and p(°°) used in (|3,100|) exist, since the 
determinants of these matrices are 1. For p(°°) see Remark 13 .17| for P n see (|3.75|) and Remarks 
l3~T£l and l3~2Tl and finally, for P n see (13371) and Remarks ET231 and I3~2%1 

By definition, R has jumps on the contour However, in the next proposition we will 
show that R has only jumps on the reduced contour see Figure 



£"' U £3 U E4" U dU 5 u at/a. 



(3.101) 



Proposition 3.31 The matrix valued function R defined by 1)3.100(1 is analytic in C\E^, where 
Y,r is defined by (|3,101j) . see a/so Figure® 

Proof. By construction, the jumps of S and p(°°) agree on Eg', the jumps of S and P ra agree 
on Xyj=i^') an( i the jumps of S 1 and P n agree on Ylj=i^j- Therefore, R has an analytic 
continuation to C \ (E^ U {0, 1}). It now suffices to show that the isolated singularities of R at 
and 1 are removable. 
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Since the Airy function is bounded near 0, it follows from (|3.59j) that P n is bounded near 
1, and since detP n = 1 thus also P^ 1 - Furthermore, S remains also bounded near 1. So, by 
(|3.100[) . R has a removable singularity at 1. 

It remains to prove that the possible isolated singularity of R at is removable. We work as 
in |191 I21| using the behavior of R near 0. This will be determined by multiplying the behavior 
of S near (which follows from (|3.52|) and condition (d) of the RH problem for T) with the 
behavior of P" 1 near (which follows from (|3.97[) together with the behavior of ^> a near 
given by j211 equations (6.19)-(6.21)] and the fact that detP n = 1). After a straightforward 
calculation we obtain, 



R(z) 
R(z) 



O 



q | log z log z 
log 2 z log 2 z J ' 



as z — > 0, if a < 0, 

as z — ► 0, if a = 0, 



(3.102) 
(3.103) 



and, if a > 0, that 



R(z) = < 



O 



o 



1 1 



as z — » 0, from outside the lens, 
, as z — > 0, from inside the lens. 



(3.104) 



From this it follows that, in all cases, R has a removable singularity at 0. If a < 0, this is clear 
from (|3.102|) and ()3.1U3|) . since a > —1. If a > 0, it follows from (|3. 104(1 that R remains bounded 
if we approach from outside the lens. Therefore, R cannot have a pole at 0. Furthermore, we 
also get from (|3.104j) that z k R(z) is bounded near for any integer k > a. Then, R cannot 
have an essential singularity at either, so that is a removable singularity of R. This ends the 
proof of the proposition. □ 

From (|3. 100(1 and condition (c) of the RH problems for S and p(°°), one then obtains that 
R is a solution of the following RH problem on the contour S^. 



RH problem for R 

(a) R : C \ E R -> C 2x2 is analytic. 
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(b) R+(z) = R-(z)v R (z), for z G E R , with 
' P n (z)P^(z)~\ for z G 8U S , 

M z ) = \ P n (z)P(°°){z)- x , for z G d&s, (3.105) 

P(°°)(z)t;5(z)p( 00 )(z)- 1 , for z G Ef U Eg' U Ef , 
where is the jump matrix of S given by (|3.53[) . 

(c) R(z) = I + 0(1/ z), as z -> 00. 

Now, we will briefly explain that R defined by ()3.100|) is the unique solution of the RH 
problem for R and that R is uniformly close to the identity matrix as n — > 00. Again we follow 
HE! Section 7.2]. 

Introduce the matrix valued function Ar = vr — I on T,r. This matrix satisfies the following 
estimates as n — > 00 (with c\ > some constant and E w = E"' U S3" U E4'), cf. jlOl Proposition 
7.7], 

\\^R\\iMfiU s udUs) + \\^R\\L 2 (dU s udU s ) + W^R\\L°°(dU s UdU s ) = C '( 1 / n )' (3.106) 
II^kIIl^e'") + II^-rIIl 2 (E"') + ||ArIU°°(£"') = C(e~ nci ), (3.107) 

uniformly for 5 in compact subsets of (0, <5o) and for v £ ^). Here, (j3.106j) follows from 
equations (|3~7l)|) and (|33S|) . Estimate (|3.107j) follows from equations (|515|l . (|3~1B1) and (JXMJ)- 

Let C_ be the Cauchy operator on E# given by C_/ = (Cf)-. The estimates for A# above 
imply that the integral operator Ca r defined as 

CaM = C-C/Afl), for / G L 2 (Er, C 2x2 ), 

is a bounded linear operator from L 2 (Y>r, C 2x2 ) into itself with operator norm HCa^ || = 0(l/n), 
as n -> 00. Therefore, Id — Ca r can be inverted by a Neumann series for n sufficiently large, 
and we define, 

fiR = (H-Ca^-HC-Ar) g L 2 (E*,C 2x2 ). (3.108) 

As in |10l Theorem 7.8] one can then show that for n sufficiently large the RH problem for R 
has a unique solution, and that 

R = I + C(Ar + ^ r Ar). (3.109) 

Note that Ar is exponentially small on E //; . Therefore, the contribution to hr from E'" 
is exponentially small. Furthermore, from (|3.7fij) . (|3.98j) and the facts that £ n and (fi n can be 
expanded in powers of n" 1 /" 1 , it follows that Ar possesses on the disks an aymptotic expansion 
in powers of n~ 1 ^ m . This will imply that /j,r possesses an asymptotic expansion in powers of 
n — i/m_ gi m ii ar as j n jim Theorem 7.10], this discussion then leads to the following theorem. 

Theorem 3.32 The matrix valued function R has the following asymptotic expansion in powers 

Of ri~m } 

^ 00 

R(z) ~ / H — y rk(z)n~™, as n — > 00, (3.110) 
k=0 

uniformly for 5 in compact subsets of (0, 5q), for v G ^-) and for z G C \ E^. Furthermore, 
the scalar functions r& are bounded functions which are analytic in C \ (dUs U c^L^) and which 
can be computed explicitly. 

This theorem states that R is uniformly close to the identity matrix as n — > 00. By going 
back in the series of transformations Y 1—* U \—* T S 1— > R we then find the asymptotics of Y. 
This ends the asymptotic analysis of the RH problem for Y. 
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4 Asymptotics of the recurrence coefficients a n ,&„_i and the 
leading coefficient 7 n 

In order to determine the asymptotics (as n — > oo) of a n , b n -\ and j n we will make use of the 
following result, see for example Let Y be the unique solution of the RH problem for Y. 

There exist 2x2 constant matrices Y\ and Y2 such that 

and 



^=(% + |r, fen-l = \ZfnR%, 7n = W 9 ■ (4.2) 

(*i)i2 V -2vrz(yiji2 

We will now rewrite the above expressions for a n , b n _\ and j n in terms of the solution of the 
rescaled RH problem for U. From (|3,14|) and (14,1(1 we have, 



_ (/ + « + * + Oil/*)) a 8 2 - co. 

Therefore, U has an expansion of the form (|4.1|) at infinity. Further, with U\ and U2 the 
analogues of Y\ and Y2, respectively, we have, 



Yl — PnPn U\p n , — P n Pn ^2Pn 



Inserting this into (|4.2|) we arrive at, 

a n = Pn (Vl)ll + T^f) . b n-i = /W(^i)i2 (^1)21, (4.3) 



and 



We thus need to determine the constant matrices U\ and XJi. For large \z\ we have by (|3.37|) . 
(13321) and (|3.1U»[) . 

So, in order to get U\ and C/2 we need the asymptotics of p(°°\z), e™ 9 "' 2 '' 73 and R(z) as z — > 00. 
Asymptotics of P(°°)(z) as 2 — > 00: 

From (|3.54|) and (|3.55|) it is easy to check that that the Szego function D has an asymptotic 
expansion in powers of z^ 1 given by D(z) = 2~ a (l + \olz~ 1 + ^a(a + 3)z~ 2 + 0(z~ 3 )), as 
z — > 00. Clearly, the scalar function a, defined by (|3.58j) . has also an asymptotic expansion in 
powers of z _1 given by a(z) = 1 — \z~ x — ^z" 2 + C(z" 3 ), as z — > 00. Inserting these asymptotics 
into (l3~5?l) we find after a straightforward calculation, 

z 

with 



p(°°)( z )=I + ^— + ^— + 0(l/z 3 ), asz^oo, (4.6) 



/ _a i\ / a(q-3)+l ■ (a+2) \ 

p(°°) _ 2 _acr 3 I 4 4 ] 2 QlT 3 p(°°) _ 2 _a<7 3 I 32 16 j 2^0-3 y\ 

1 " I _i « / ' 2 ~ I -(a-2) a(a+3)+l I ' ^ ' ' 

V 4 4/ \ «^6" 32 / 
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Asymptotics of e ngn ^ a ' A as z — > oo 
By (|3.31j) we have, cf. equations ( 

ng n (z)a 3 ( 2 



and 





z 11 



with 



G^ 



-n Jq ti> n (t)dt 







n Jq tip n (t)dty 



Ch, 



as z — > oo, 



* 
* 



(4- 



(4.9) 



Asymptotics of -R(-z) as z — > oo: 

Analogous as in ^U] the matrix valued function R has the following asymptotic expansion at 
infinity, 



R(z) =I+?± + ^ + 0{l/z 3 ), asz^oo, 
z z z 

where the 2x2 constant matrices R± and i?2 satisfy, 



(4.10) 



Ri 



Ri 



1 

1 

2ni 



dU 5 



(y)dy 



dU s 



yAi(y)dy 



1 

2tH 



2-ki 



A l {y)dy)- + 0{l/n I 
du s I n 



dU s 



1 



yA 1 (y)dy) - + 0(l/n 



n 



as n — > oo, 



as n — > oo, 



with Ai and Ai given by (|3.7T|) and (|3.99j) . respectively. Recall that Ai is meromorphic in 
{0 < \z — 1| < 5o} with pole of order 2 at 1, and that Ai is meromorphic in {0 < \z\ < So} with 
pole of order 1 at 0. So, one can evaluate the above integrals using the residue theorem. After 
a rather long (but straightforward) calculation we obtain, 



Ri 



-aa s 



1 



+ 



16M0) 

1 



(4a 2 



1) 



1 



1 



+ 



h'M) 



16^(1)2 



1 -i 
-i -1 



-12a 2 + 3 i(12a 2 + 24a + 11) 
48/i n (l) VjK 12 " 2 - 24 « + 11) I 2 " 2 - 3 



and 



12 



i4~ 



12a 2 + 24a + 16 _ h' n (l) \ 1 2a 
48^(1) ^ ' ~ +U[l/n 



16/i n (l) 2 / n 



2 a<T3 - + C(l/n 2 ), (4.11) 
n 



(4.12) 



Remark 4.1 The error terms in (|4.11|) and (|4.12l) have an asymptotic expansion in powers of 
n -V m w ith coefficients that can be explicitly determined, cf. |10| Section 8.1]. 



We now have all the ingredients to determine the asymptotics of a n , b n -\ and j n . Note that 
by (|4.5|1 . (|4.6j) . (|4.8() and (|4.1U|) the constant matrices U\ and U2 are given by 



Z7i = 



(4.13) 



and 



U 2 = e ^ nln(TZ 



P M + G2 + p 2 + + +R 1 )G 1 



e -\nt n oz 



(4.14) 
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Proof of Theorem 12. 1L First, we determine the asymptotics of b n -\. Inserting (j4.13j) into 
(14,31) . and using the fact that (Gi)i2 = (^1)21 = 0, see (|4,9[) . the recurrence coefficient 6 n _i is 
given by, 



p (oo) )i2(p (oo) )2i + (Pi ^ )l2(jRl ) 21 + (R l)l2 (P^>) 21 + (i? 1 ) 12 (i? 1 ) 21 



,(oo) 



(00) N 



1/2 



From equations (|4.7I) and (|4.11|) . and from the fact that (^1)12(^2)21 = 0(l/n 2 ) as n — > 00, we 
then arrive at 



%-i 



A, 

4 

4 



[1 + 4* (4~ a ( J R 1 ) 2 i - 4 Q (i?!) 12 ) + 0(l/n 2 



4a 1 ,„ . , 9 N 
1 + OITn +0 ( 1 /" ! ' 



1/2 



as n — > 00. 



(4.15) 



This proves equation l|2.5|) . 

Next, we consider a n . Inserting (j4.13|) and 1)4. 14|) into (|4.3|) . and using the facts that (Gi)i 2 
(C 2 )i2 = and (Cn)ii + ((^1)22 = 0, see (j4.9|) . the recurrence coefficient a n is given by, 



a n = Pn[ (Pi (oo) )n + (#1)11 + 



(P2 (00) )l2 + (^2)l2 + (^lP 1 (00) )l2 , 

>(°°)> 



\ (pr')i2 + (Ri)i2 j 

Now, from ()4.7j) and from the fact that (i?i)i2 = 0(l/n) as n — > 00, it follows that 
1 



(4-16) 



(i'rtia + (#1)12 



-4i4 a (1 + 4i4 a (i?i)i2 + C(l/n 2 )) , as n -> 00. 



Inserting this into (j4,16|) . we then find, from (|4.7j) . (|4.11j) and (|4.12|l . after a straightforward 
calculation, 



ftu 

2 



a n = ^(l + 4(i2i)n + 4i4 Q (i?!) 12 - 8i4 a (i? 2 )i 2 + 0(l/n 2 )) 
2(a + 1) 1 



1 + 



/t n (l) n 



+ 0{l/n I 



as n — > 00, 



which proves (|2.6|) . 

Finally, we consider the leading coefficients j n . Inserting (|4,13|) into (|4.4() . and using ()4.7j) . 
(|4.9|) and the fact that (i?i)i 2 = 0(l/n) as n -> 00, we find, 



7n = A 



-(n+f+l)„-in«n./ 2 



e 2' 



-2 Q (l-4i4 Q ( J R 1 ) 12 ) 

7T 



-1/2 



I3 n {n+ i + i ] e -ht n ^_2 a (1 + 2*4 a (fli)i2 + 0(l/n 2 )) , as n -> 00. 



(4.17) 



From (|4~TT|) we then obtain (|2~7)) . 

Note that the statement that all the error terms have an explicit asymptotic expansion in 
powers of n _1 / m follows from Remark 14.11 This ends the proof of Theorem 12.11 □ 



5 Plancherel-Rotach type asymptotics for the orthonormal poly- 
nomials p n 

In this section we will determine the asymptotic behavior (as n — * 00) of p n (Pn,z) in the four 
asymptotic regions Ag, B$, C$ and D$, see Figure ^ We will do this by rewriting p n {(3 n z) in 
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terms of the solution U of the rescaled RH problem for U. From (J32J, (jSHH), (|3~T2j) and (j22J), 

we obtain 

Pn(Pnz) = (Jte)-%e&^e-^rJ^ ^e-^y^Uniz)) (1 + 0(1/0), (5-1) 

as n —* oo, uniformly for z G C. 

So, we need to determine the asymptotics of U\\. In order to do this we will make use of the 
following properties. Since g n — £ n — \i n — ^V n is analytic in C \ (— oo, 1] and zero on (1, oo), 
see (|3.44[) . it follows from the uniqueness theorem that 

9n{z)-£n{z)-\tn = \Vn{z), for Z G C \ (-OO, 1]. (5.2) 

The reader can easily verify that z l l 2 + [z — l) 1 / 2 = (p(z) 1 / 2 and z 1 ! 2 — (z — 1) 1//2 = ^p(z)^ 1 ^ 2 
for z G C \ (— oo, 1], where 99 is defined by (|2.11l) . From (|3.58f) we then obtain, 

a(z)+a(z)~ 1 _ V {z) 1 / 2 a(z)-a(z)- 1 _ . if{z)-^ 2 



2 " 2zV*(z - l) 1 / 4 ' _2i 2zV4(z - 1)1/4 ' 

for z G C \ [0, 1] by analytic continuation. Inserting this into i|3.57j) and using the fact that 
<p(z) = e i™<™>s(2z-l) for z G C+, we have 

2*V4(2 - 1)1/4 y_. ^ z) \{a-X) ^ z j-I(o-l) J 



z -2 aCT 3 ; for z G C+, (5.3) 



2z x l A (z - I) 1 / 4 I _j e M»(*) e -*»eW 
where we have introduced (for the sake of brevity) the notation, 

771(2:) = -(a + 1) arccos(2z — 1), V2(z) = —(a — 1) arccos(2z — 1). 
5.1 The outside region As 

In this subsection we determine the asymptotics of the orthonormal polynomials p n in the outside 
region A$. 

Proof of Theorem E3] (a). For z G A s we have from (|3~37|) . (|3~52|) . I|3.100|> . from the first 
equality in (|5.3|) . and from ()5.2I) . 

[/(*) = . ^ „ . r-f e ^w e ^iz(z)2- ag8 f v^) |( ; +1 M (54) 

W V / 2zV 4 (z - I) 1 / 4 V; \-iip{z)^ a - 1 ) ) V 7 

So, from Theorem 13.321 and from the fact that 1/c/? is uniformly bounded in As, we then obtain 

U n (z) = zie^^We^^-^^e^W (R n (z) - i4 a ' R l2 {z) V {z)- 1 ) 

= z -^ e hnV n {z) e \ni n2 -« ^00 2(Q+ ^ e ng„( Z ) ^ + 0(1/ n )) , 

2z i ' 4 (2: — l) 1 / 4 

as 71 — > co, uniformly for 5 in compact subsets of (0, 5$) and for z E As. Inserting this into 1)5. 1|) 
and using the definition of £ n , see (|3.41[) . the first part of the theorem is proven. □ 
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5.2 The bulk region Bs 

Here, we will determine the asymptotics of p n in the bulk region Bs- We will make use of the 
following proposition, which will also be used (in the next section) to prove the universality 
result in the bulk of the spectrum, cf. Theorem 12.81 fa). 

Proposition 5.1 For z E Bs, the first column of U satisfies, 

U{z) Q = z-f e 1 2 nV "^e 1 2 ne ^R(z)P^\z)z c i^ (J^^j , (5-5) 

where £ n is given by (|3,41|) . p(°°) is the parametrix for the outside region given by 1)3. 57|) . and 
R is the solution of the RH problem for R given by 1)3.100)1 . 

Proof. The proposition can easily be verified using 1)3.37)1 . 1)3.52)1 . 1)3.100)1 and 1)5.2(1 . □ 

Proof of Theorem 12.41 (b). Inserting the expression 1)5.3(1 for into equation 1)5.5(1 . and 

using the fact that (2— l) 1 / 4 = (1 — z) l / A e~ for z E Bs, and the fact that £ n (z) = —^i J* ip n (s)ds, 
we arrive at, 

U{z) (^]= __l__ 2 -f e |nK( Z ) e |n^3 



oy zV4(i_ z )i/ 



(cos f 77i (-2) — ff ib n (s)ds — t) \ 
wn > Ji vn\ > 4 ; (5 6) 
-icos (r? 2 (^) - vrn / x ip n (s)ds - \ )J 



The second part of the theorem now follows easily from this equation together with 1)5.1(1 and 
Theorem 13.321 □ 

5.3 The Airy region C$ 

In order to determine the asymptotics of the orthonormal polynomials in the Airy region Cs 
we start, as in the previous subsection, with a proposition which will also be used to prove the 
universality result at the soft edge of the spectrum, cf. Theorem 12.81 (b). 

Proposition 5.2 For z E Us, the first column of U satisfies, 



U(z) 



J) = yfoe-Zz-U^We^WEnW (a^S)) ' (5 ' 7) 

where f n is the scalar function given by ()3.65() . E n is the 2x2 matrix valued function given by 
1)3.74(1 . and R is the solution of the RH problem for R given by ()3, 100(1 . 

Proof. Let z E Us be such that f n (z) E II, cf. Figure El This means that z lies inside the disk 
U s and belongs to the upper part of the lens. Then, from (EOTl) . (13321 . (13.1001) . ((3~75l and ()PJl 
we find after a straightforward calculation, 

£/(*) Q = z-f e^ v '"( z ) e l^ CT3 J R(z) J B n (z)M/(/„(z)) Q , (5.8) 

where the 2x2 matrix valued function ^ is given by ()3.63() . Now, since f n {z) E II, we have by 
(I3~fi3l . 

Inserting this into ()5.8|) . equation ()5.7() is proven in this sector of Us- The proof in the other 
sectors (i.e. z E Us such that / n (z) E I U III U IV) is analogous. □ 
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Proof of Theorem 12.41 (c). Inserting expression 

Q for P(°°) into the definition (|3~74l of 

E n , we find 



z 1 / 4 (z-l) 1 / 4 v^—i cos 772(2) -sin 772(2 

for z£ Inserting this in turn into (|5.7|) we arrive at, 



U{z) I J) = V^r -f e inV»( z ) e |n^ 3 



0; zV4( z _ 1)1/4' 

w \-icosr] 2 (z) -smr] 2 (z)J Jny ' \Ai' (f n (z)) J ' v 7 

for z G C<5. Now, one can easily check that the scalar functions cos 772/ cos 771 and sin 772 / sin 7/1 
are bounded in Cg, so that by Theorem 13.321 

Uu(z) = z -" le hV n {z) e \ni n2 - a _ 



7T 



31/4(3-1)1/4 

x (cos7 ?1 (z)(l + 0(l/n)) -iBtotfrOOa + Oa/n)))/^ (5.10) 

as n — ► 00, uniformly for S in compact subsets of (0, So) and for z G C5. Inserting this into (|5,1|) 
part (c) of the theorem is proven. □ 

5.4 The Bessel region Dg 

We will start with a result similar to Propositions I5.ll and 15, 2l This result will also be used to 
prove the universality result at the hard edge of the spectrum, cf. Theorem 12. 81 (c) and Theorem 

Proposition 5.3 For z G Ug, the first column of U satisfies, 

U(z)(^ = z-fe^e^R(z)E n (z) ( ^'t^Tl . ^ J , (Ml) 



-2vri(-/„(z))i/ 2 J4(2(-/ n (z))V2) 



where f n is the scalar function given by (|3.91[) . i£ n is i/ie 2x2 matrix valued function given by 
(|3.96j) . and i? is i/ie solution of the RH problem for R given by ()3.100|) . 

Proof. Let z G % be such that f n (z) G III', cf. Figure H This means that z lies inside the disk 
U s and belongs to the upper part of the lens. Then, from (EP71) . (|3~52l . (jH.lOQI) . (13371) . (t5~2l . 

and from the fact that (— z)~2 ao " 3 = z - 2 aaz e^ ama3 for z G C+, we find after a straightforward 
calculation, 

U(z) (J\ = z-^e^ nVn ^e^ nl ^R{z)E n {z)^ a (/„(z)) e^ maa ' A (J\ , (5.12) 

where the 2x2 matrix valued function *$> a is given by ([3.82ft . Now, since f n (z) G III', we have 
by (|3~%2l and Q] formulae 9.1.3 and 9.1.4] 



*a (/„(*)) e5™ CT 3 A) 



J Q (2(-/ n (z))V2) 
2 7 r/ n (z) 1 /V4(2(-/ n (z)) 1 /2) 



Inserting this into (|5.12|) and using the fact that f n {z) 1/2 = -i(f n {z)) 1/2 for z G C+, the 
proposition is proven in this sector of Ug. The proof in the other sectors (i.e. z G Ug such that 
f n (z) G Full') is analogous. For the case f n {z) G I' one has to use formulae 9.1.35 and 
9.6.3]. □ 
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Proof of Theorem 12.41 (d). Inserting (j5,3j) into the definition (|3.9fij) of E n , and using the 
facts that for z G D$, 

= z^e-^, (z - I) 1 / 4 = (1 - z) 1 ^, Uz) 1 ' 4 = (-/ n (z))V4 e -f ) 



we obtain 

Ejz) 



-l) n 2" 



00-3 



sinCi(.) ico8Ci(*)w ; n ( z ))^( 2vr )¥ 



y/2z l / A (l _ z )l/4 ^-isil^C*) C0SC2(^) 



(5.13) 



for z G Ds- Here the scalar functions Ci and C2 are given by Ci,2(^) = ^1,2(2) — vtq/2. Plugging 
this in into 115.111). we arrive at, 



K ) \0J { iJ Z V4(l-^)l/4 

x R{z)2~ aa3 



Smdiz) COsCl(^) \ (Ja{2{-fn{z)) l/2 i 

-isinQ 2 {z) -icosQ 2 {z)J 1 j/ f ( 2 (-/„(z)) 1 / 2 ) 



, (5.14) 



for z £ Ds- Now, since the functions sin £2/ sin£i and cos £2 / cosCi are bounded in Ds, we then 
have by Theorem 13.321 



_a 1 

C/n(z) = z 2 e 2 



ny„(z) |n«„9-a/ -iw V / ^(-/n(^)) 1/4 
6 2 Zl/4(1 -z)l/4 



/ 4 (1-Z)1 

x (sinCi(^)(l + 0(l/n)) cosCi(^)(l + 0(l/n))) 



/ J Q (2(-/ n (z))V2)^ 
J' a (2(-f n (z))^) / 



, (5.15) 



as n — > 00, uniformly for 5 in compact subsets of (0, <5o) and for z G D$. The last part of the 
theorem then follows from (|5.1|) . □ 

Proposition 15.31 gives the behavior of the first column of U near the origin. To get the 
universal behavior of the three kernels Wi >n , Wii jTl and Wm> at the the hard edge, see Theorem 
I2.1U1 we also need the behavior of the second column of U near the origin, cf. Remark l3.3l This 
will be given by the next proposition. 

Proposition 5.4 For z G Us, the second column ofU satisfies, 



U(z) 



z^e 2 



\nV n (z) hni n a 3 



e* n ^ a3 R(z)E n (z) 



\H^(2{-f n (z)) 



1/2 > 



x < 



-7ri(-/ n (z)) 1 / 2 (F«)'(2(-/ n (z)) 1 / 2 ) 

-i^ 2) (2(-/ n (z))V 2 ) 
ni(-f n (z)y/ 2 (H^)'(2(-f n (z)) 1/2 ) 1 
Proof. The proof is analogous to the proof of Proposition 15.31 



ifze u s nc + , 



if z G u s n C_ . 



(5.16) 



□ 



Remark 5.5 In work in progress |SJ we need to evaluate U(0) for the case a > 0. From the 
fact that -f n {z) = c n n 2 z(l + 0(z)) as z ->• 0, and from (l5~TTIl . (t5~T6l . (ETT^ and 1, formulae 
9.1.7, 9.1.9 and 9.1.27], we obtain, 



U{0) = e2 ne ^R(0)E n (0) 



1 

"2vr 



Cn n a 

W) 



us 



C2 



5<9(0)<T3 
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So, it remains to evaluate E n (0). This can be done using (|5.13j) and the facts sin Ci (z) = l + 0{z), 
sinC 2 (-z) = -l + 0(z), cosCi(-z) = (a + l)y/z(l + 0{z)) and cosC2(-z) = (1 - ct)yfz(\ + 0{z)) as 
z — > 0. After a straightforward calculation we then arrive at, 



17(0) = {-l) n e^ nlnaz R{U)2- aa * (\ * j**^) 











(i 













OS 

I • (5.17) 



Note that det f7(0) = 1 as it should be. 

6 Applications in random matrix theory 

In this final section we will prove Theorems 12.81 and I2.1UI We will only consider the proof of 
Theorem 12.81 (c). i.e. the universality of K n at the hard edge of the spectrum, in detail. The 
other results can be proven analogously, and are left as an easy exercise for the reader. 

We start by writing the kernel K n in terms of the first column of U. Recall that K n is 
already written in terms of the first column of Y by (|3.5j) . Using (|3.14j) we then obtain, 

n K n (f3 n x,(3 n y) = x% e-f^-Wyf e -^(v)_J_ ( l) ir\y)U(x) Q . (6.1) 

Then, we prove Theorem 12. 81 (c) by inserting the expression for the first column of U derived in 
Proposition 15.31 into this expression for K n . In order to prove part (a) and (b) of the theorem 
we have to use Propositions 15 . ll and 15 . 2l respectively. 

Introduce, for the sake of brevity, some notation. With c n = (|/i n (0)) 2 , cf. Remark 13.261 
define 



u v 



■/„W) 1/2 , and v n = 2{-} n (v ri )) 1 ' 2 . 



Ac n n 2 '' n 4c n n 2 ' 
Note that by fl3~32]) and (l3~Ml) . 

u n = \/u(l + 0(u/n)) and v n = y/v(l + 0(v/n)), as n — > oo, (6-2) 

uniformly for u,v in bounded subsets of (0, oo). Further, we also need the following lemma. 

Lemma 6.1 Let L n = E n R, where E n is given by ()3.96j) . and where R is the solution of the 
RH problem for R. Then, with the notation of above, 

L~ 1 (u n )L n ('u n ) = I + - , as n ► oo, (6.3) 

uniformly for u, v in bounded subsets of (0, oo). 

Proof. Let u, v € (0, oo) and let 7 be a closed contour in Us encircling the origin once in the 
positive direction. Then, since L n is analytic in U§ we have by Cauchy's formula, 

L-\v n )L n (u n ) =I+±-.l L n 1 K)^n(C) _ ^ 



2ttz J 7 (C -u n )(( 

for u, v in bounded subsets of (0, 00) and n sufficiently large. Now, from (|3.96ft and from the 
fact that R is uniformly (in n and z) bounded in Ug, we have that L n {Q) = 0(n 1//2 ), as n — > 00, 
uniformly for £ G 7. Furthermore, since detL n = 1 we obtain that L~ 1 (w n ) = ©(n 1 / 2 ), as 
n — ► 00, uniformly for v in bounded subsets of (0, 00). Together with u n — v n = O(^), this 
proves the lemma. □ 
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Now, we have introduced the necessary ingredients to prove the universal behavior of the 
kernel K n at the hard edge of the spectrum. 

Proof of Theorem 12.81 (c). With the notation of above, we have by (|5.11j) . 

U(u n ) = u^e^ v ^e^ e ^L n (u n ) ( J ?^}_ ,) ■ (6.4) 
\0J \-mu n J' a (u n ) J 

Furthermore, since A f = (\ J) A' 1 (5V) for any 2 x 2 matrix A with det A = 1, we obtain 
(note that det U = detL n = 1), 



(0 l)U-\v n ) 



v J e hnV n (v n ) frivnJ'^) J a (v n ))L- l {v n )e-^\ (6.5) 



Inserting ()6,4j) and 1)6.5)) into (|6.1|1 . we find from 1)6.3 J) and from the facts that J a (u n ) = 
0{u a / 2 ) = u n J' a (u n ) and J a (v n ) = C(W 2 ) = v n J' a (v n ) as n -> oo, 

= 55b) (i » J « (5 » ) J " (i "» (-5;!)) + (^) ■ (6 - 6) 

as n — > oo, uniformly for u, -u in bounded subsets of (0, oo). We can now replace u n and v n by 
-y/u and y/v, respectively, see ()6.2)1 . and we then arrive at the Bessel kernel. However, then we 
make an error which could not be estimated uniformly for u — v close to zero. So, we need to be 
more careful, and we refer the interested reader to the presentation in 22 , Proof of Theorem 1.1 
(c)] where this problem was solved (note that ()6.6|) corresponds to [HI equation (3.31)]). Then, 
part (c) of the theorem is proven. □ 

Proof of Theorem 12.81 (a) and (b). The proof of the first two parts of the theorem is 
analogous to the proof of part (c). For part (a) we have to use ()5.5)) . and for part (b) we use 
()5.7)) . Further, the scaling with c n n 2 / 3 at the soft edge of the spectrum is clear from ()3.7U)) . □ 

Proof of Theorem 12.101 The proof is analogous to the proof of Theorem 12.81 We first have 
to write (like we have done in ()6.1j) with K n ) the kernels W\^ n , W\x >n and Wni,n in terms of U. 
From H3,6j) and (|3.14|) we have 

PnWu tn (P n V,f3 n u) (3nW U l,n(PnU, (3 n v) 
-PnWl tn ((3 n U, (3 n v) -(3nWu tn (f3 n U, f3 n v] 

~ 2^7— -pt^U-'^U^P'^ 3 . (6.7) 



Then, we prove the theorem by inserting the expressions for the first and second column of U 
(derived in Propositions 15.31 and IB~4l respectively) into this expression. This is left as an easy 
exercise for the reader. □ 

Remark 6.2 Since Wi >n depends only on the first column of U, it is clear (by Proposition 15.31) 
that Wi >n depends only on the J-Bessel functions. Further, since W\i^ n depends on both the first 
and second column of U, this kernel depends (by Propositions 15.31 and 15.4)) on both the J-Bessel 
and Hankel functions. And finally, since VFn^n depends only on the second column of U, it is 
clear (by Proposition 15.4)1 that this kernel depends only on the Hankel functions. 
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